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_^ Abstract 

o 

O \ We present a detailed description of a zero temperature phase transition between 

;_i ■ superconducting and diffusive metallic states in very thin wires due to a Cooper 

^r<, pair breaking mechanism. The dissipative critical theory contains current reducing 

c/j \ fluctuations in the guise of both quantum and thermally activated phase slips. A 

full cross-over phase diagram is computed via an expansion in the inverse number of 
complex components of the superconducting order parameter (one in the physical 
case). The fluctuation corrections to the electrical [a) and thermal (k) conductivities 
'^ ' are determined, and we find that a has a non-monotonic temperature dependence 

in the metallic phase which may be consistent with recent experimental results 
CJ , on ultra-narrow wires. In the quantum critical regime, the ratio of the thermal 

to electrical conductivity displays a linear temperature dependence and thus the 
Wiedemann-Franz law is obeyed, with a new universal experimentally verifiable 
J^ . Lorenz number. 

I> 
(N 

O. 

2? ■ 1 Introduction 

> 

r> ' At the nanoscale, the basic mechanical, electrical and optical properties of 

c^ ' materials that are well understood at macroscopic length scales can change 

in interesting and sometimes unexpected ways as quantization and fluctua- 
tion effects manifest themselves. An intriguing question thus arises regarding 
the implications of reducing the scale or effective dimensionality of materi- 
als, that even in the bulk, are known to already display interesting quantum 
mechanical behavior. Superconductors are natural physical systems to con- 
sider in this context. Conventional or low temperature superconductors are 
well understood in the bulk, unlike their high temperature cousins whose full 
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description still remains elusive after more than twenty years of intensive re- 
search. A major obstacle to the study of high temperature superconducting 
materials is that they are plagued by their proximity to competing states with 
both order and disorder at the atomic scale. Through a better understanding 
of the ways in which normal superconductivity is suppressed or destroyed in 
different confining geometries and effective dimensions, perhaps progress can 
be made towards a mastery of this fascinating emergent phenomena at all 
length and temperature scales. 

Recent technological advances now allow experimentalists to fabricate and 
control systems with ever smaller dimensionality. In two dimensions, thin films 
of lead with less than ten atomic layers can be grown with such precision that 
mesas and valleys are intentionally engineered to trap vortices, hardening the 
superconductor in a magnetic field ([l[). In one dimension, a technique known 
as suspended molecular templating (}2|) can produce long metallic wires with 
diameters less than 10 nm. Finally, the fabrication and ultimate measurement 
of spin excitations in manganese chains using a scanning tunneling micro- 
scope in inelastic tunneling mode is an elegant example of the observation 
and control of a system consisting of only a few atoms (J3|) . 

In this paper, we will focus on the transition between a superconductor and 
a metal in ultra-narrow wires. It is well known that the Mermin-Wagner- 
Hohenberg theorem (|j; |5|) precludes the possibility of long range supercon- 
ducting order at any non-zero temperature in one dimension. However, any 
real wire is three dimensional, and can be approximated as a cylinder with a 
finite radius R. In the Landauer picture (jGl) conduction is proportional to the 
number of channels in the wire, A^_l, equal to the number of states that can be 
occupied for a given energy in all dimensions transverse to transport. Thus, if 
we imagine free electron states propagating down the wire, A''_l oc A/\p where 
A is the cross-sectional area and Ap is the Fermi wavelength. Any real wire will 
have R^ Xf implying that A^_l ^ 1 and we can imagine that it will undergo 
a phase transition to a superconducting state below some critical tempera- 
ture. As the diameter of the wire decreases, or at suitably low temperatures, 
it will eventually enter a regime where the superconducting coherence length, 
^ equal to the average separation between Cooper pairs, is larger than the 
radius. This condition defines the quasi-one dimensional limit, as paired elec- 
trons necessarily experience the finiteness of the transverse dimension while 
unpaired electrons do not. 



1.1 LAMH theory 



The study of superconducting fluctuations in narrow wires has a long history 
beginning in 1968, when Webb and Warburton performed a remarkable trans- 



port experiment on thin whisker-crystals of Sn (|7|) with diameters between 40 
and 400 /im. They noticed that for the thinnest wires, resistive fluctuations 
leading to a finite voltage persisted below the bulk critical temperature (T^) 
for tin and defied any mean field characterization using the maximum super- 
current (jsl). The understanding of this behavior, which is specific to quasi-one 
dimensional superconductors followed rapidly thereafter and is composed of 
three parts. The first was Little's qualitative introduction of thermally acti- 
vated phase slips (J9|); an unwinding of the phase of the superconducting order 
parameter by ±27r in a region of the wire where the magnitude of the su- 
perconducting order parameter has been spontaneously suppressed to zero. 
These non-trivial thermal fluctuations are equivalent to a vortex tunneling 
across the wire. In an applied bias current, negative jumps in the phase can 
exactly balance the linear in time phase increase needed by the Josephson re- 
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lation for the system to exhibit a finite voltage below T^ (llOl ). Next came the 
Ginzburg-Landau (GL) theory of Langer and Ambegaokar (jlll ) for the free 
energy barrier height of a phase slip event which qualitatively reproduced the 
most important features of Webb and Warburton's experiments. The story 
concluded with the time-dependent GL theory of McCumber and Halperin 



(Il2l ) who correctly computed the rate at which these resistive fluctuations 
occur, and led to full quantitative agreement. The contributions of Langer, 
Ambegaokar, McCumber and Halperin are now referred to as the LAMH the- 
ory of thermally activated phase slips. 

Transport near the finite temperature phase boundary between the supercon- 
ducting and normal state is controlled by these phase slips, which from the 
LAMH theory have a free energy barrier height AF{T) and occur at a rate 
0(T) at temperature T. When used in conjunction with the Josephson rela- 
tion, LAMH showed that their contribution to the resistance of the wire is 
given by 

^^(^) ^-AF(T)/fcBr n T\ 

where Rg = h/{2ef' is the quantum of resistance. This result can be applied 
to a nanowire of length L and normal resistance R^ = {4:L/N_j_i)Rg where Nj_ 
is the number of transverse channels and i is the mean free path to obtain the 



LAMH contribution to the conductivity (J13l ) 
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or 



e^ fT N±i\ 

C^LAMH = -J^^^^ *LAMH ( ^, 77^ 1 (1-3) 

where $lamh is a universal dimensionless function, T^ is the superconducting 
transition temperature and ^(0) is the zero temperature GL coherence length. 



A multitude of experiments on superconducting nanowires (J2|; [ij; [l3|; [l5|; [l6 



171 : Il8l : Il9l : 1201 ) have confirmed the accuracy of the LAMH theory by fitting 



Eq. (I1.2p to experimental transport measurements with T^ and ^(0) as free 
parameters with great success. This description includes the effects of only 
thermally activated phase slips, and neglects the possibility of quantum phase 
slips at low temperatures (l2ll:l22l) where, if present, one would expect quantum 
tunneling to produce deviations from the LAMH theory. There are some exper- 
imental indications that quantura phase slips (QPT) may indeed be present at 
the lowest temperatures (1231 : 124 : |25| ) and we attempt to address some of these 
issues in Section [23] by presenting a version of the LAMH theory with param- 
eters renormalized by quantum fiuctuations. For a recent and comprehensive 
review describing the infiuence of both thermally activated and quantum phase 
slips on superconductivity in one dimension see Ref. (1261 ). 



1.2 Ultra-narrow wires 



As the diameter of a wire is reduced, there are two important changes that 
need to be considered. The first is the well known volume to surface area 
ratio, and thus surface effects will begin to affect bulk behavior. The second is 
more subtle and is related to the increased effects of coupling to an external 
environment. In the presence of such dissipation, a small system can undergo 
a quantum localization transition, as is observed in small Josephson junctions 

(0). 



In the early 1990s, step edge electron beam lithography techniques were used 
to create narrow indium strips with diameters between 40 and 100 nm. When 
transport measurements were performed, there appeared to be significant devi- 
ations from the LAMH resistance at low temperatures resulting in a persisting 
resistance manifest as a "foot" raised upwards from the expected exponentially 
decreasing resistance (J23l : |24J ). It was proposed that this was due to the onset 
of quantum phase slips at low temperatures occurring via the macroscopic 
tunneling mechanism of Caldeira and Leggett (J28|). These ideas were vigor- 



ously pursued (l29l : l30l : l31l : |32| ) leading to a host of theories which did not 



necessarily agree on the observability of quantum phase slips in experiments. 
One of the most interesting results was an upper bound on the wire diameter 
of approximately 10 nm above which quantum phase slips would be strongly 
suppressed ( i3(^ ) as their rate ^qps ~ exp(— iV^) can be exponentially small, 
where N±_ is the number of transverse channels in the wire discussed above. 
This upper bound of ten nanometers was far too narrow for step edge lithog- 
raphy techniques and it would take the invention of new fabrication methods 
to fully address these issues. 



Wires with truly nanoscale dimensions were not studied until the introduc- 



tion of a novel and pioneering nanofabrication technique known as suspended 
molecular templating in early 2000 (l^). This remarkable process can be used 
to manufacture wires with lengths between 100 and 200 nm with diameters 
less than 10 nm. The key feature is the top down approach that uses a long 
narrow molecule such as a carbon nanotube or DNA as a backbone on top of 
which the wire is deposited. The fabrication process begins by etching a trench 
in a substrate formed from a silicon wafer using electron beam lithography. 
The backbone molecules are then placed in solution and deposited over the 
substrate. They are allowed to settle, and at high concentrations some will end 
up resting over the trench. The entire surface of the substrate is then sputter 
coated with several nanometers of a metal like Nb or alloy such as MoGe. The 
result is that a thin uniform layer of the deposited material is suspended over 
the trench by the backbone molecule. It can be located via scanning electron 
microscopy (SEM) and then isolated with a mask that is also used to pattern 
electrodes that will be used for transport measurements. 

When measuring the resistance of a given wire as a function of tempera- 
ture, two exponential dips are observed. The first at high temperatures, cor- 
responds to the large two dimensional leads going superconducting while the 
lower temperature drop is due to the actual wire undergoing a transition. The 
temperature at which the wire goes superconducting is strongly dependent on 
its diameter, with thinner wires being pushed to lower temperatures. Superb 
agreement with the LAMH theory is found for normal state resistances down 
to 0.5 Q using Eq. (11. ip . At resistances below this value, or for thinner wires, 
there appears to be a growing experimental consensus that there is a quali- 
tative change in the resistance including deviations from the theory of purely 
thermally activated phase slips. The wires seem to be entering a regime where 
resistive fluctuations coming from other effects, possibly including Coulomb 
blockade and quantum phase slips can either postpone or completely destroy 
the superconducting transition (2|; [l3|; [l9|; l33l ) . This behavior can be seen by 



measuring the resistance of thinner and thinner wires as a function of temper- 
ature leading to a separation between superconducting and metallic transport 
all the way down to the lowest temperatures as seen in Fig. [H If supercon- 
ductivity is indeed being destroyed as the temperature is reduced to zero by 
quantum and not thermal fluctuations upon tuning some parameter related 
to the size of the transverse dimension, then such a transition is by definition 
a superconductor- metal quantum phase transition (1351 ) (SMT). 



1.3 Pairbreaking quantum phase transition 



From the Cooper instability in BCS theory it is known that the existence of 
a non-trivial quantum critical point in a metal implies a finite electron inter- 
action strength. Naively the existence of interactions would seem to preclude 
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Fig. 1. Experimental transport measurements on MoGe nanowires reproduced from 
Ref. (J34l ) showing a distinct difference between thick superconducting and thin 
metalhc or resistive wires as the temperature is reduced. At zero temperature, a 
quantum critical point would separate the superconducting and metallic phase, with 
the transition between them being described by a quantum superconductor-metal 
transition (SMT). 



the possibility of any quantum phase transition between a superconductor 
and a metal because the ordered phase will always exist at some strictly non- 
zero temperature for arbitrarily weak paring. Moreover, if the temperature is 
driven to zero in a pure BCS superconductor, pair fluctuations will be com- 
pletely eliminated (J36l ) . The solution arrives in the form of pairbreaking inter- 
actions, or any perturbation that is odd under time reversal symmetry. Such a 
term will act differently on the spin and momentum reversed constituents of a 
Cooper pair, making pairing more difficult. The presence of these interactions 
effectively cuts off the logarithmic singularity in the pair susceptibility and 
sets a critical value for the strength of the pair potential before superconduc- 
tivity can develop. Therefore, our proposed SMT must live in the pairbreaking 
universality class. 



The mean-field theory for the SMT goes back to the early work (1371 ) of 
Abrikosov and Gor'kov (AG). In one of the preliminary discussions of a quan- 
tum phase transition, they showed that a large enough concentration of mag- 
netic impurities could induce a SMT at T = (T^ is protected by Anderson's 



theorem (138l ) for the case of non-magnetic impurities). The transition is tuned 
by a parameter a which is proportional to the impurity concentration and AG 
derived an equation for the phase boundary given by 

where '\\){x) is the polygamma function, k^T^a = l.lAhiUDe"^^^^'^^^ is the BCS 
transition temperature in the absence of any pairbreaking perturbations with 
ud the Debye frequency. Eq. (11. 4p shows that by perturbing a conventional 
superconductor with a suitably strong interaction that breaks time reversal 
symmetry, it is possible to completely destroy the superconducting state at 
finite temperature at adT). Mathematically, this is equivalent to the obser- 
vation that for large enough a, Eq. (11.41) has no non-zero solution for the 
temperature. 

It has since been shown that such a theory applies in a large variety of sit- 
uations with pairbreaking perturbations fi39i ). For the general case, ha can 
be interpreted as the depairing energy or splitting between the two time- 
reversed electrons of a Cooper pair, averaged over the time required to com- 
pletely uncorrelated their phases. Relevant examples include anisotropic gap 



superconductors with non-magnetic impurities (J40l:l4ll:l42l). lower-dimensional 



superconductors with magnetic fields oriented in a direction parallel to the 



Cooper pair motion (J43l: |44|) . and s-wave superconductors with inhomogeneity 
in the strength of the attractive BCS interaction (45). Indeed, it is expected 
that pairbreaking is present in any experimentally realizable SMT at T = 0. 
In the nanowire experiments, explicit evidence for pairbreaking magnetic mo- 
ments on the wire surface was presented recently by Rogachev et al. (20). Any 
impurity at the surface would be much less effectively screened, and one could 
imagine a BCS coupling V{r) which depends on the radial coordinate of the 
wire r. In this picture, V{r) would change sign from negative (attractive) to 
positive (repulsive) as r changes from r = to r = i? where R is the diameter 
of the wire. This behavior is schematically outlined in Fig. [2l For the thickest 
wires, R> C, and the mean field solution to the BCS equations will lead to the 
wire being described by a superconducting core surrounded by a cylindrical 
metallic envelope. A similar picture was recently put forward to describe the 
SMT in two dimensions by imagining superconducting grains embedded into 
a film with a pairin g in teraction which depended on the distance from the 



center of the islands (|45| ). In the nanowire case, as the transverse dimension is 
reduced, there will be a critical radius, R ^ C,, where the superconducting core 
will vanish and the wire will enter a metallic state. This is a rather physically 
appealing picture as it is suitable for the destruction of superconductivity in a 
wire that is only weakly disordered in the bulk and is well suited to theoretical 
models. 
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Fig. 2. A schematic cross-section of a metallic wire where magnetic impurities on 
the surface are poorly screened leading to a change in sign of the BCS pairing 
interaction as one moves from the center, |^| 7^ to the edge, |^| ~ where ^(r) 
is the superconducting order parameter. Below some temperature, the wire would 
be composed of a superconducting core with a normal resistive sheath. 

It is clearly impossible to continually reduce the diameter of a single wire while 
measuring its transport properties and thus other more systematic examples 
of a SMT would be beneficial. The most obvious candidate is a transition 
that can be observed in a single wire by increasing the strength of a magnetic 
field oriented parallel to its long axis. Such an experiment has been performed 



by Rogachev, Bollinger and Bezryadin (1l7l ) on individual Nb nanowires. As 
the strength of the parallel magnetic field is increased, the superconducting 
transition appears at lower and lower temperatures with the expectation that 
for suitably strong fields it will vanish all together and the wire will exhibit 
metallic behavior; a quantum SMT. To completely destroy superconductivity, 
pairbreaking events must be uncorrelated over long time scales and a theoret- 
ical description would require diffusive electrons or suitably strong boundary 
scattering. 



1.4 Microscopic approach 



Fluctuations about the AG theory have been considered (l43l : |44J ) in the metal- 
lic state, and lead to the well-known Aslamazov-Larkin (AL) (l46l) . Maki- 
Thompson (MT) (1471 : l48l ) and Density of States (DoS) (|49| ) corrections to 
the conductivity. The form of these corrections is usually introduced in terms 
of the structure of their diagrammatic representation within the finite tem- 
perature disordered electron perturbation theory but they all have the same 
physical origin: in the presence of strong pairbreaking, the normal metal still 
experiences pairing fluctuations near the Fermi surface as a result of its prox- 
imity to the superconducting state. Specifically, the AL effect comes from the 
direct charge transfer from fluctuating Cooper pairs, the MT correction results 
from coherent Andreev scattering off the fluctuating pairs and the Density of 
States (DoS) correction is due to the reduction of the normal electron density 
of states near the Fermi surface accounting for the paired electrons. 



The exact form of these contributions are known from recent microscopic 
finite temperature perturbative computations in BCS theory (J43l : |44| ). These 
results are vahd at low temperatures, with the pairbreaking parameter a. larger 
than critical Uc of the SMT. Denoting the BCS coherence length by ^o? one 
can define a clean limit by ,^0 *^ ^ and a dirty limit by £ <^ ,^o- The total 
conductivity was obtained in the dirty limit and it was found that 
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where ctq is a background metallic conductivity, c is a non-universal con- 
stant, D is the diffusion constant in the metal, and the remaining correc- 
tions from pairing fluctuations have been written in the form of a power of 
T times a factor within the square brackets which depends only upon the 
ratio h{a — ac)/kBT. Writing the conductivity in this way allows us to de- 
termine the relative importance of the fluctuations corrections, in the renor- 
malization group sense, to the SMT. The first square bracket represents the 
usual Aslamazov-Larkin (AL) correction (J46l ) and has a prefactor of a nega- 
tive power of T and so is a relevant perturbation. This is a result of the large 
inverse power oi a — ac which will be dominant near criticality as the critical 
point is approached. The second term arises from the additional AL, MT and 
DoS corrections: the prefactor has no divergence as a power of T, and so this 
correction is formally irrelevant at the SMT. The complete second term has 
a finite limit as T — > 0, and so becomes larger than the formally relevant AL 
term at sufficiently low T in the metal. The second term is therefore identified 
as being dangerously irrelevant in critical phenomena parlance: it is important 
for the properties of the low T metallic region, but can be safely ignored at 
finite temperatures near the critical coupling. 



1.5 Field-theoretic approach 



At the SMT, field-theoretic analyses (J50l : l5ll ) show that the AG theory, along 
with the AL, MT and DoS corrections, is inadequate in spatial dimension 
d < 2, and additional repulsive self-interactions among Cooper pairs have to 
be included. Here, d defines the dimensionality of the Cooper pair motion. The 
confining dimension, or radius of the wire, R, is larger than the inverse Fermi 
wavevector, but smaller than the superconducting coherence length or Cooper 
pair size, ^. This is the exact condition discussed previously for the quasi-one 
dimensional limit. While the Cooper pairs are effectively one dimensional, any 
unpaired electrons have a three dimensional Fermi surface and thus strictly 
Id Luttinger liquid physics do not apply. 
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Fig. 3. Crossover phase diagram of the superconductor-metal transition in a 
quasi-one dimensional wire. The metal (M) is described by the perturbative theory 
of Ref. (|43l ). The quantum critical (QC) region is described by our effective theory 
Sa for temperatures above Tdis where the effects of disorder may be neglected. The 
Mooij-Schon mode is present everywhere, but couples strongly to superconduct- 
ing fluctuations only in the fluctuating superconductor (FSC) regime, where it is 
described by Eq. ()2.24|) : note that Sa does not apply here. The dashed lines are 
crossover boundaries which occur at T ~ |a — ad (from Eq. (j2.8p ) and the interme- 
diate LAMH region should be described by the theory of thermally activated phase 
slips discussed in Section 11.11 with the possible modifications of Chapter [3l Note 
that Sa does not include the MT and DoS corrections in the metallic regime where 
they are actually larger than the AL contribution to the conductivity. 

The goal of this work is to understand the aforementioned experiments on 
nanoscale metallic wires and we are interested in the fluctuation corrections to 
the thermal and electrical conductivity across the SMT as well as the nature 
of the crossovers from this universal quantum critical physics to previously 
studied regimes at low T about the superconducting and metallic phases. In 
the remaining sections, we will examine the d = 1 SMT in great detail and 
begin by introducing a theory for an ultra-narrow superconducting wire in 
terms of a strongly-coupled field theory of bosonic Cooper pairs, overdamped 
by their coupling to fermionic quasiparticles (normal electrons). Both the zero 
frequency thermal and electrical transport coefficients are computed, first in 
the large- iV limit, where the number of components of the superconducting 
order parameter (A^) is assumed to be infinite, and then as an expansion in 
1/A^. Along the way we make contact with the microscopic BCS theory, touch 
upon an extension of the LAMH theory near T^ and investigate the ratio of 
the thermal to electrical conductivity, known as the Wiedemann-Franz ratio. 

To frame our discussion, we conclude the introduction with a summary of 
the full crossover phase diagram that is constructed throughout this paper 
in Fig. [3l The important features of this global pairbreaking phase diagram 
include a strongly fiuctuating quantum critical regime, the main focus of this 
study, where phase and amplitude fiuctuations must be treated on equal foot- 
ing. For a pairbreaking strength greater than the critical one, as the tempera- 



10 



ture is reduced, there is a crossover to a low-T metallic regime described by the 
theory (J43l : |44J ) of AL+MT+DoS corrections in d = 1 (dashed line in Fig. [3]). 
On the superconducting side, there is a regime of intermediate temperatures 
where the classical LAMH phase slip theory applies (I 111 : Il2l ). with possible 
modifications coming from quantum renormalized coefficients and eventually 
another crossover at still lower temperatures to a phase fluctuating regime 
whose description requires a non-linear a-model of fermion pair fluctuations 
coupled to the superconducting order (1521) . 



Models similar to the one we analyze, but lacking amplitude fluctuations (so 
called phase-only theories) have been previously applied to the physical sit- 
uation considered here (J30l : l53l : |5J; |55|; l56l ) . In these continuum theories, the 
destruction of superconductivity is due to the proliferation of quantum phase 
slips resulting in a normal phase which we maintain is insulating and not 
metallic at T = 0. Within the phase only approach, the superfluid's conductiv- 
ity would be controlled by irrelevant phase-slip operators as the resulting d = 1 
superconductor-insulator transition (SIT) lives in the Kosterlitz-Thouless uni- 
versality class (ISTI ). Such theories would thus be applicable to a quantum SIT 
which could be appropriate in short inhomogeneous wires. 



2 Dissipative model 



The app roach to the SMT taken here is akin to previously studied theories 



(J58l : l59l : l60l ) of disordered superconducting fllms with unconventional pair- 



ing symmetry. Such fllms are assumed to be composed of a network of small 
Josephson coupled superconducting islands. The transition between the nor- 
mal and superconducting state can be tuned by altering the distance between 
the grains; if the separation becomes large enough, quantum fluctuations can 
destroy superconducting order even at zero temperature. Physically, the tran- 
sition occurs when the Josephson coupling becomes less than the Coulomb 
energy due to the transfer of a single Cooper pair between islands. These ar- 
guments can be made more rigorous by starting from microscopic BCS theory 
and deriving an effective model for the fluctuations of the Cooper pair order 
parameter. The order parameter is strongly overdamped by decay into quasi- 
particle excitations manifest as an interaction that is long range in imaginary 
time (0; M, H H H S) • 



These ideas can be directly applied to the d = 1 quantum superconductor- 
metal transition in the pairbreaking universality class which is relevant for 
ultra-narrow quasi one-dimensional metallic wires. The physical ingredients 
include repulsive paired states and a lack of charge conservation of the con- 
densate requiring the existence of a bath (the large number of transverse 
conduction channels) into which the former member electrons of a disassoci- 
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ated Cooper pair can flow into upon breaking. The final description will be in 
terms of a strongly-coupled field theory of bosonic Cooper pairs, overdamped 
by their coupling to the unpaired fermionic states of the metal. 

The fluctuation corrections to transport associated with the AL correction 



discussed above are naturally captured in this picture. These have (1431 : l44j ) 
a Cooper pair propagator {Dk^ + \ujn\ + a)~^ at wavevector k and imagi- 
nary bosonic Matsubara frequency iun in the metal in both the clean and 
dirty limits. Here, the "mass" or bare pairbreaking frequency a measures the 
strength of the pairbreaking interaction which could come from a variety of 
sources as mentioned in the introduction. D is equal to the usual diffusion 
constant D = D = vpi/3 [vp is the Fermi velocity) in the dirty limit where 
the mean free path is much smaller than the superconducting coherence length 
{£ ^ ^o)- In the clean limit, where i ^ C,o, D will be in general some non- 
universal number that depends on the specific microscopic details (such as 
the lattice constant) of the system in question. This motivates the quantum 



critical theory of Ref. (J50l : l6ll ) for a field \E'(x,r) which represents the local 
Cooper pair operator 



L l/T 

Sa = dx dr 





2 , „,I,T,^^ ,M2 , ^|,Tr/™ ^M4 



D\d.,^{x,T)\' + a\^{x,T)\' + -\^{X,T 



1^ 

+ TY,jdx'y\uJn\\'^{x,uJn)[' 



(2.1) 



where we have used the temporal Fourier transform of \E'(x, r) 

l/T 

■^{x,ujn) = f dr '^{x,T)e''^"\ (2.2) 



with uJn = 2TTnT to more compactly express the non-locality of the dissipative 
term in imaginary time and have chosen units where h = ks = 1 for conve- 
nience. From this point forward, we will suppress the limits of integration for 
the sake of compactness unless their inclusion is required for clarity. The quar- 
tic coupling u must be positive to ensure stability, and describes the repulsion 
between Cooper pairs. The pairs are strongly overdamped, and the rate of 
their decay into the metallic bath is characterized by the coupling constant 
multiplying the \un\ term, 7, which is required to be positive by causality. It 
will be convenient to rescale the field \E' such that the coefficient of the Landau 
damping term is equal to unity. In addition, we rescale all couplings according 
to 

^^ ; D ^ 'jD; a ^ 'ja; u -^ 'j'^u. (2.3) 

This theory describes the vicinity of a superconductor-metal quantum critical 
point, corresponding to the (bare) value z = 2 for the dynamic critical expo- 
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nent. A different description {z ^ 2) cannot be completely ruled out, but it 
would most likely require additional tuning parameters as well as the inclusion 



of unusual pairing phenomena (1361 ) . 



The quantum phase transition is driven by altering the strength of the pair- 
breaking frequency a as was shown schematically in Fig. [31 While a ^ ac 
there is normal metallic conduction and for a -^ ac the system is fully super- 
conducting. For a > ac pairing fluctuations enhance the conductivity while 
for a < Oc, both thermal and quantum phase shps, included as amplitude 
fluctuations of \E' that destroy the superflow. 

The field theor y in Eq. (12.1 p is identical in form to the Hertz-Millis-Moriya 



theory (J62l : l63l : l64l ) describing the Fermi liquid to spin-density wave (SDW) 
transition, with the Cooper pair operator \E' replaced by an 0(3) order param- 
eter representing diffusive paramagnons. In the neighborhood of this transi- 
tion, k measures the magnitude of the deviation from the SDW ordering wave 
vector K and the dissipative |ct;„| term arises from the damping of order pa- 
rameter fluctuations resulting from coupling to gapless fermionic excitations 
of the metal near points on the Fermi surface connected by K. A more careful 
analysis leads to the realization that at T = on the ordered (SDW) side of 
the transition, a gap appears in the fermion spectrum for small k. Thus, this 
description is only fully accurate at T = on the disordered (metallic) side of 
the transition or at finite temperature in the quantum critical regime above 
the SDW state. The same logic applies to the role of phase fluctuations at low 
temperatures in the superconducting phase near the SMT. We will return to 
this point later with a thorough discussion of the Mooij-Schon normal mode, 
but for now we begin with a detailed scaling analysis of Sa- 



2.1 Scaling analysis 



Simple power counting for the rescaled self-interaction term in Eq. (12.11) in d 
spatial dimensions fixes the upper critical dimension at c? = 2 and necessitates 
a non-perturbative treatment for the quasi-one dimensional limit considered 
here. In d = 1, m is a relevant perturbation but in the strong coupling regime 
{u -^ oo) we expect all results to be universal {u independent). 



Pankov et al. (l65l ) studied a theory similar to S^ with \1' replaced by a A^- 
component field via the renormalization group (RG) in an e = 2 — d expansion 
in one and two dimensions at zero temperature. The most important result 
obtained from their RG analysis is that the damping term, \uJn\-i generated from 
the long-range 1/t^ interaction between order parameter fluctuations does not 
require an independent renormalization, and thus the frequency dependence 
of the propagator only involves wavefunction renormalization. At T = and 
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a = Oc in one dimension, they find a non-trivial fixed point and an analysis 
of the RG equations leads to an expression for the dynamical susceptibility at 
small frequencies and momenta 



X{k, u)= JdrJ dxi^lix, r)^,(0, 0))e- 



i(kx—ujT) 



where a = 1, . . . ,N, $^ q is a universal scaling function and Cq is a non- 
universal constant that will depend on D. The dynamical critical exponent 
can thus be read off as 

z = 2-r] (2.5) 

where its bare value has been corrected by an anomalous dimension rj ~ 
e^. This result holds to all orders due to the existence of only wavefunction 



renormalization mot . l66l ). Eq. (12.41) can be generalized to finite temperatures 
where 

x(A:,-„,T) = i$,(^,|l^), (2.6) 

with $^ another universal scaling function and Ci a non- universal constant. 
Most interestingly, at fc = and ujn = the value of the inverse susceptibility 
in the quantum critical region will be fixed by temperature alone, 

X-\0,0) = AT (2.7) 

and the highly non-trivial universal constant A will be computed in a 1/A^ 
expansion in Section HI 

Scaling functions for the most singular parts of the dc electrical a and thermal 
K, conductivities can also be derived with a knowledge of their scaling dimen- 
sions alone. In one dimension, the longitudinal conductivity of a wire is equal 
to e^/h times a length, and for /c = a;„ = but finite temperature there is 
only one length scale available, the thermal length Lt ~ T~^/^, and the energy 
scale is set by the distance from the critical point. Similar arguments apply 
for the thermal conductivity leading to (returning to physical units) 




where z/ is the usual correlation length exponent defined by ^ ~ |a — ad"'^ 
and $o- and $« are two dimensionless universal scaling functions. 

In the Gaussian (non-interacting) limit, z = 2, rj = and u = 1/2 and 
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corrections can be computed in the e = 2 — d expansion (J65l : ISOl ) 

(A^ + 2)(12-7r 



4(Ar + 8)2 



-e' + 0{e') (2.10) 



_ 1 (A^ + 2) (iV + 2) [6A^2 ^ (228 - 77i'^)N + 792 - SStt^] ^ 3 

'^ ~ 2 ^ 4(iV + 8)^ ^ 48(Ar + 8)3 ^ ^ ^^ ''" 

(2.11) 

These results are in agreement with Monte Carlo simulations (67. ) which found 
z = 1.97(3), z + T]= 1.985(20) and u = 0.689(6). 

If interactions are included but calculations are performed in the limit where 
the number of order parameter components is large, the large-N limit, z and 
7] are unchanged from their Gaussian values, but v = 1 and the replacement 

[n{a - fte)]^ M« - «c) . ^ 

is required in the scaling functions above. Corrections from the N = 00 values 
of the exponents u and 77 can be computed in the 1/A^ expansion and the 
results are detailed in Section HI 



The RG analysis of Ref. (J65l ) confirms that Sa satisfies conventional hyper- 
scaling relations at the T = SMT in the absence of disorder. This implies 
that all irrelevant operators can be neglected, and transport should be fully 
described by Eq. (12. ip . Moreover, the most singular part of the dc conductiv- 
ity of a given wire will be described by Eq. (12.81) which is independent of its 
length L. 



2.2 Particle-hole asymmetry 



In the scaling analysis of the previous section, we neglected an important 
detail; for z = 2, when the energy dependence of the electronic density of 
states near the Fermi level is included, a propagating term in the action can 
arise from the weak particle-hole asymmetry of the electronic spectrum 

Sp= I dx I dTp^>*{x,T)—m{x,T). (2.13) 

The magnitude of particle-hole symmetry breaking, p is proportional to the 
energy derivative of the density of states at the Fermi energy, and we therefore 
expect its bare value to be small. This is supported by the fact that our quasi- 
one dimensional treatment of Cooper pairs coupled to a bath composed of 
three dimensional electrons, required the ratio of the pairing to Fermi energy 



15 



to be small. Specifically, if Ap ^ -R ^ ^, then using the BCS result for the zero 
temperature energy gap A = e-p I i^-K^k-p) we find that l/fcp^ ~ A/ep ^ 1. 



In fact, a microscopic weak coupling derivation of Eg. (12.1) fo r a dirty two 



dimensional superconductor with d-wave pairing symmetry (l40l : l42l ) finds that 
the ratio of the dissipative to propagative terms is proportional to the dimen- 
sionless conductivity of the normal phase, p oc l/epT where r is the scattering 
time in the self-consistent Born approximation. For a good metal, the product 
of the scattering time and Fermi energy is large, and thus p <^ 1. This is con- 
sistent with results at finite temperature for a weak-coupling short coherence 
length superconductor (JGSl ) . 

At tree level, p is marginal, and as just argued, we expect its bare value to be 
small. However, we can examine the renormalization group fate of Sp near the 
fixed point of Sa- The scaling dimension of p can be computed m. a. d = 2 — t 
expansion for a massless (a = 0) quantum critical theory through the conven- 
tional method of isolating any logarithmic singularities (or 1/e poles) in the 
Feynman diagrams corresponding to all possible insertions of the perturbing 
term, ipoj in Eq. fl2.13p . For the case considered here, there are two unique 
graphs 



X \ (2.14) 





where a solid line is equal to the bare propagator (fc^ + \uj\) ^, a dot represents 
the quartic interaction u and a cross is an insertion coming from Eq. ( 12.13p . 

If the external lines have frequency VL and zero momentum, then the combi- 
nation of these graphs leads to the integral 

'^""^ = -'^P"" J ^ J ^ J (2^J 



{27rrJ (2vr)^ (2.15) 

"" (F + |a;i|)2(g2 + \uj,\)[(k + qf + \uj^ + CU2 + fi|] ' 

A simple power counting analysis of the integrand in d = 2 — e dimensions 
leads to the appearance of the predicted pole 

/(fi)=zf]p^, (2.16) 

where as usual, the flow equation for p is related to the residue A via 

%-A,. (2.1T) 
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I{Q) is computed in Appendix |X] leading to the flow equation 

dp V? 



IGvr^ V TT 



1-- p. (2.18) 



The fixed point value of u is given in Ref. (J6S) for the equivalent z = 2 0{N) 
model with the change of notation uq = 3u. Equivalently it can be easily 
computed to one loop order as 

2-71^ 

u* = -^e (2.19) 



fvr^ 



V)p (2.20) 



leading to 

dp ^ 

~dl ~ 100' 

or p[£) ~ e^''^^ ^ at RG scale I. Thus we conclude that although p is relevant, 
its scaling dimension is extremely small. In conjunction with a bare value that 
we have argued should be diminutive, we will neglect Sp in future calculations. 

There is still one piece missing in our analysis of S^ as alluded to in the 
previous section; the role of charge conservation (which S^ breaks) and the 
associated normal modes. 



2.3 Phase fluctuations 



From hydrodynamic arguments, it is known that a one-dimensional metal or 
superconductor should support a gapless plasmon, or a Mooij-Schon normal 



mode (1691 ) ■ which disperses as u; ~ k\n^''^{l/{kR)). Our discussion of this issue 



parallels that in Refs. (J70l : l7ll ) on the role of conservation laws in the critical 
fiuctuations of quantum transitions in metallic systems for which the order 
parameter is overdamped (as is the case here). To make this explicit, couple 
\E' to a fiuctuating scalar potential A^- with bare action 

^ rdk r dujAr{k^uj)f_ 
^ J 27tJ 27r A\n{l /{kR))' ^ " ^ 

However, the nature of the A^--^^ coupling differs between the "quantum criti- 
cal" and "fluctuating superconductor" regimes of Fig. [31 For the main results 
of this study, we need only the coupling in the quantum critical region, where 
the physics of the plasmon mode is unchanged from that in the "Metal" region 
of Fig. [3l After integrating out the fermions, we obtain the A^ action 

Su = J^J^lll{k,u^)\Mk,u.)\', (2.22) 

where 11 is the irreducible density correlation function (the "polarizability" ) 
of the metal. For u ':^ k, we have Il{k,Lu) ~ k'^/uj'^, and then Sa + Sn has a 
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pole at the plasinon frequency noted above. We also observe that the coupling 
between these A^ fluctuations and \1/ is negligible: a slowly varying A^. is a shift 
in the local chemical potential, and to the extent we can ignore the variation 
in the electronic density of states at these energy scales, the effective couplings 
in Sa do not change with Ar, and there is no Ar — "^ coupling. The existence 
of an Ar — "^ coupling that is non-analytic in frequency has been discussed by 
loffe and Millis (70'). However, they show by Ward identities, which also apply 
here, that these couplings do not contribute to the physical charge correlations. 
So, in the quantum critical region, Sa and Sa + Su are independent theories 
describing the pairing and charge fluctuations respectively. 

As described in the previous section, the coupling constant which measures 
particle- hole asymmetry is formally, albeit weakly relevant (see Eq. fl2.20p ) and 
thus a small A,- — \1/ coupling appears by making the r-derivative in Eq. (12.131) 
Gauge covariant 



Sp = I dxdr 



P'^*{x,t) l—-2eiAr] ^ix,T) 



(2.23) 



However, the combination of the small scaling dimension and the small bare 
value of 7] implies that such particle-hole asymmetric effects, and the con- 
sequent coupling between pairing and charge fluctuations, can justifiably be 
ignored in theories hoping to describe realistic experiments. 

We conclude by addressing the physics in the "Fluctuating superconductor" 
regime of Fig. [3] for a < a^- In this discussion, we neglect the possibility of 
a narrow region a ^ a^ with gapless superconductivity, which is known to 
occur in the vicinity of a pair breaking transition (1391) . Now coupling between 
the pairing and charge fluctuations is much stronger. When T < (oc — a) the 
action S^ does not apply for the smallest frequencies. The reasons for this are 
again analogous to arguments made for the spin- density- wave ordering tran- 



sition in metals, as discussed in Section [2] and Ref. (J35l ). For the latter case, 
it was argued that with the emergence of long-range spin density wave order, 
the low energy fermionic particle-hole excitations at the ordering wavevector 
were gapped out, and so the diffusive paramagnon action applied only for 
energies larger than this gap. At energies smaller than the gap, spin- waves 
with dispersion u; ~ fc emerge. In the superconducting case, there is no true 
long-range order at any T > 0, but the order is disrupted primarily by 'renor- 
malized classical' thermal fluctuations of the phase, </> of the complex \1/ field. 
We assume that there is a local pairing amplitude in the fermion spectrum, 
analogous to the spin-density wave order. The low energy effective action for 
obtained by integrating the fermions in the presence of a local pairing, is 

S^ = JdxdT{[Ki [dr(l)ix,T) - 2eArf + K2 [d,(l)ix,r)f} (2.24) 

where K12 vanish as power of {ac — a) (.71. ). The strongly coupled pairing 



and charge fluctuations in the "Fluctuating superconductor" regime of Fig. [3] 
are described by Sa + Sfj,, and this theory contains the Mooij-Schon mode, 
which is the analog of the 'spin-wave' mode. We do not claim that Sa + S^ 
can extended across quantum criticality into the normal phase, in contrast to 
other works (30. ) which consider vortex unbinding in such a theory. 

The arguments presented in this section strive to justify the use of Sa through- 
out the quantum critical regime where the presence of repulsive interactions 
between Cooper pairs will be essential in obtaining the correct form of the full 
crossover transport scaling functions $o- and $k- This is most clearly seen by 
examining the energy and length scales at which various terms in Sa become 
dominant in both the clean limit where the mean free path is longer than the 
superconducting coherence length, and the dirty limit, where the opposite is 
true. 



2.4 Universality and interactions in the quantum critical regime 



The microscopic theory of the superconductor-metal transition was considered 
in great detail by Shah and Lopatin (SL) (44i) in a Gaussian theory of super- 
conducting fluctuations that corresponds to the effective fleld theory presented 
here with u set to zero. In Section [3] we will compare the transport proper- 
ties computed from the effective action Sa with Shah's results in the metallic 
regime where a ^ ttc- Before doing so we will need to connect the (until now) 
phenomenological coupling constants of Eq. (12. ip to those of microscopic BCS 
theory. The details are given in Appendix [B] with the main result being that 
in the dirty (d) limit, the diffusion (D), dissipation (7) and interaction (u) 
constants are given by 

Dd = D = \v^i (2.25) 

u, . 2.9^ (2.27) 

where the number of transverse conduction channels in the wire A^^ is assumed 
to be large and in the clean (c) limit 

Dc = \vFio (2.28) 

2.0 

u, = Ud- 2.9-^. (2.30) 
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Note that the microscopic value of the quartic couphng constant u is identical 
in both the clean and dirty limit. 

Armed with the microscopic values of the model parameters we can make a 
number of observations regarding the validity of the non-interacting theory 
in the strongly fluctuating quantum critical regime (which is referred to as 
the classical regime in SL). In the non-interacting theory, the temperature 
dependence of the superconductor-metal phase boundary adT) is computed 
from an expansion of the mean-field result of Abrikosov and Gor'kov given in 
Eq. (11.41) for T <^ a^- In SL, the deviation from criticality is measured with 
respect to the finite temperature mean field phase boundary 

ac{T)=ac- , (2.31) 

3nTco 

where 7 ^ 0.577 is the Euler-Mascheroni constant and Tco is the classical 
BCS transition temperature in the absence of pairbreaking. This notation 
(which we temporarily adopt) requires comment. The location of the quantum 
critical point is as usual defined as Oc at T = 0, but in Eq. (12.311) a^T) is 
the function which locates the value of the pair breaking frequency at which 
superconducting order is lost for a given temperature. It is the approximate 
functional inverse of Eq. (II. 4p in the low temperature limit. To summarize, SL 
are interested in large positive values of a far into the metallic phase and have 
chosen to define a coupling constant that measures the distance from classical 
criticality defined by the temperature dependent mean field phase boundary 
and not the distance from quantum criticality. 

To successfully compare the approach considered here with that of SL, we 
shift the definition of a accordingly and write (returning to physical units) 



dx / dr 
kuT 



D\d.^i^,r)\^ + ^T'\^ix,r)\' + ^Mx,r)\' 



H ^Xl / ^^l^"ll^(^'^") 



(2.32) 



In this form, it is clear that the theory only goes quantum critical at T = 0. The 
coupling constants D and u can take on the values computed in Appendix [B] 
for the clean and dirty limits, but we are primarily concerned with the role 
of the quartic coupling u, characterizing the strength of the Cooper pair self 
interaction. For d = 1, u has scaling dimension one and from Eq. (IB.lSp it 
has engineering dimensions of inverse mass times inverse length or frequency 
squared times length over energy. There are three distinct regions of the phase 
diagram as the temperature is reduced in the quantum critical regime set by 
the size of the bare value of u defined by the conditions: 

I The quartic coupling can be ignored and the Gaussian theory of SL fully 



20 



describes transport. 
II Interactions are important, and Hartree corrections to the mass must be 

included leading to non-universal results. 
Ill The quartic coupling is relevant, its bare value is large, and all results 
are universal. 



Cases I and II can be distinguished by examining the lowest order correction 
to the perturbatively renormalized or Hartree corrected mass coming from 
Eq. (I2.32P at one loop order 



R 



ShTco 

dk 



duj 



1 



27rDF+|u;| 



h 



T.T 



Dk"^ + \uJn\ + Ti'ykBT'^/ShTcoJ 

(2.33) 



where we have applied the usual shift to subtract off a zero temperature con- 
tribution so that our renormalized mass i? = at quantum criticality. The 
Hartree correction can be separated into two contributions, one coming from 
the integral, and one coming from the most dominant contribution to the sum, 
the Un = term. These are given by 



du r dk 



1 



2n J 2n \Dk^ + |cu„| + 7r-fkBT^/3hTco Dk^ + \uj\ 



TT 



and 



kuT r dk 



SksT, 



cO 



TT^k^T'^ 

\ 3hDT,o 
(2.34) 



(2.35) 



h J 2nDk^ + 7i-fT^/3Tco 2\^ n-fhlD 

respectively. Provided that T < Tco, the second is the most dominant contri- 
bution, and thus the renormalized mass is significant when it is greater than 
the bare mass of Eq. (12.311) . 



hu 



SkBTco TcykBT"^ 



2 ^ njhD 3hT, 
which defines the Hartree temperature 



(2.36) 



cO 



^bTh 



'ShkBT, 



cO 



3/4 



U 



1/2 



717 



.2VD, 



(2.37) 



For temperatures above Th one can ignore the presence of a repulsive inter- 
action between the Cooper pairs, and the non-interacting results of Ref. (J44j ) 
will be accurate. 
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The temperature below which all results scale to universal values can be ob- 
tained by considering the thermal length which follows naturally from our 
previous scaling analysis Lt ~ T~^/^ or more precisely for 2 = 2 



hb 

The bare quartic coupling can be assumed to be large with respect to all 
other parameters and thus flow to infinite strength when the potential energy 
is greater than the kinetic energy, i.e., 

h^u hD 

-r>JT- ^2.39) 

This relation sets the temperature T\] below which one can safely take u -^ 00 
and obtain universal results to be 

ksTv = -^. (2.40) 

The values of the microscopic parameters computed in Appendix [B] and re- 
peated above can be used to evaluate the temperatures defined in Eqs. (12.371) 
and (12.401) which separate the regions I-II and II-III. As before, the results 
depend on whether the system is taken to be in the clean or dirty limit. 



2.4.1 Dirty Limit {^0 '> ^) 

The Hartree temperature in the dirty limit can be found by substituting 
Eqs. (EH), (IRT3D and (iRlSl) in Eq. (12:371) 



Th d = 0.83-—— ^^^^-7-, (2.41) 

(6ocA^±a)'/" 

where the single electron localization length is defined to be 

60c = N^i. (2.42) 

This temperature can be converted into a length scale, which gives a lower 
bound on lengths over which one must explicitly include Hartree corrections 

Lh,. = Om^H'ii^of^'. (2.43) 

The universal temperature scale is found from Eq. fl2.40p to be 

T„.= |^^. (2.44) 

JV_l Qoc 
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Gaussian 


no n- Gaussian 
non-universal 


non- Gaussian 
universal 


DiHy Limit (^o > ^) 


hr T ^ n '^1 ^""F 


25 hvp ^ 1 rp ^ PI Qo fep 


^B^<iiT£ 


LbI ^ 0--te„,iv^53)i/4 


JVx?,oe -'*^ - "•"•^teociVx?g)V4 


Komcli'mr/' 


0-63eioc''(^^o)'/' < i < 0.126oc 


L > 0.126oc 


Clean Limit (^q ^ ^) 


^bT > 0.96^ 


33^ < k^T < 0.96^ 


fceT < 33^ 


L < 0.51iV^/^^o 


0.51N]/^Co <l< omoN^Co 


L > 0.090iVxCo 



Table 1 

The temperature and length scales in the clean and dirty limits corresponding to 

the regions of applicability described in 1-111 for the effective action Sa- 



corresponding to length scales longer than 

Lv,d = 0.12^ioc. 



(2.45) 



2.4.2 Clean Limit (^0 < ^) 



The same analysis can be repeated using Eqs. (IB. 171) and ( IB. 191) for the clean 
limit. The Hartree temperature is given by 



Th,c = 0.96- 



hvp 



^^0 



with associated length scale 



Lh,c = 0.51iVl/'eo. 



For universal results we find 



Tt 



U,c 



33 



hvp 

NUo 



with 



(2.46) 



(2.47) 



(2.48) 



Lv,c = O.OdON^Co. (2.49) 

The results are summarized in Table [H but it is immediately clear that when 
in the clean limit, L-r > L^^c can be easily satisfied, whereas in the dirty 
limit Lt > L-[J4 would require lengths on the order of ^loc, and thus weak 
localization effects could become important. 
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We must therefore restrict our analysis to temperatures greater than those 
where disorder effects set in, which we now investigate. 



2. 5 The role of disorder 



Until now, the presence of disorder in the wire, manifest as spatially dependent 
coefficients in Sa has been neglected. This topic has already been studied in 



great detail by the authors (J72l ) and will not be focused upon here. Instead, 
we choose to consider only those temperatures above which weak localization 
effects can be safely neglected and limit the scope of our results to the non- 
random universality class. An estimate of the temperature scale T^is where 
disorder effects must be included is found by equating the thermal length 
with the localization length. This yields 

Lt = 6oc = NJ, (2.50) 

and using Eq. (IZSHD and Eq. (IRT3|) we find 



"bT.. = ^, (2-81) 

where r = d/v-p is the elastic scattering time. Tjis can therefore be made 
arbitrarily small by considering thicker or cleaner wires. 

The analysis performed in this section has provided a firm foundation for 
the applicability of the effective action Sa to the SMT in ultra-narrow wires. 
In the next section transport results are computed near this quantum phase 
transition in the limit where the number of complex components of ^ is large. 



3 Transport in the Large-A^ limit 



In order to incorporate the repulsive interaction between pairing fiuctuations, 
the effective dissipative action must be generalized from the physical case of 
describing a 1-component complex field ^ corresponding to the Cooper pair 
operator, to an A^-component complex field '^a with a = 1, . . . ,N. This section 
presents a calculation of the thermal k, and electrical a dc transport coefficients 
in the "LAMH" , "quantum critical" and "metal" regimes described in Fig. [3] 
through the application of both analytical and numerical methods. It is al- 
ways assumed, unless otherwise specified, that these fiuctuation corrections 
are the most singular terms at finite temperature resulting from the direct 
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contribution to transport due to Cooper pairs, (i.e. any subscripts on trans- 
port coefficients are suppressed). It is understood that making contact with 
real experimental measurements would require a subtraction of normal state 
transport results. 



3.1 Effective classical theory 



Generalizing the effective action of Eq. (12. ip to describe the fluctuations of a 
A^-component pairing field at finite temperature 



^a 



Sn = dx dr 



L'|9^.^a(x,r)|2 + a|^„(x,r)|2 + -|$„(x,r) 



+ TJ2 dx\uJn\\'^a{x,UJn)\'^, 



(3.1) 



where we have used the short hand notation 

|^„(x,r)|2 = |^i(x,r)|2 + --- + |^;v(x,r)p (3.2) 

\^i,^(x,T)\''^[\^,{x,r)\' + ■■■ + \^^{x,r)\f (3.3) 

for the sake of compactness. 

The aforementioned goal of calculating uniform electrical and thermal trans- 
port properties in the dc limit will require special care to ensure that the zero 
frequency limit is taken while the temperature remains finite. The universal 
properties of the SMT can be most easily accessed in the strong coupling 
regime, u ^ oo, while the ratio of m to a is held fixed. The result is a much 
simpler hard spin quadratic action which is written in Fourier space as 



T r dh 
^9 = -J2 ^\^aik,uj^)nDe + |.;„|). (3.4) 

q ^ J 27C 



along with the constraint |\l/a(a^, t)P = 1 where \l/a(fc, cUn) is the Fourier trans- 
form of \l/a(a;, r) defined by 

'^a{k,Un) = j dxjdT^a{x,T)e-'^^^-^-^\ (3.5) 

The parameter g now tunes across the quantum critical point and the quantum 
partition function is given by 

Z = jvmj)^>l5{\^>f-l)e-^^. (3.6) 



We will first investigate observables at frequencies much smaller than the 
temperature in the continuum quantum critical regime {huj <C k-BT). Here, 
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temperature plays the role of an infrared cutoff and the effects of quantum 
fluctuations can be integrated out, producing an effective classical theory with 
renormalized parameters. This can be done in the large- A^ limit by first im- 
posing the constraint |\E'a(x, r)p = 1 via a Lagrange multiplier /i and then 
integrating out all non-zero Matsubara frequencies from Z over their Gaus- 
sian action. The resulting effective action has an overall factor of A^, and as 
A^ — > cxo we can perform the functional integral over [i using the saddle point 
approximation where we replace r = i[i. This yields the classical partition 
function 

Z, = |l)vl/^I)v^;exp{-^|rfa;p|a,^,(a;)p + \/(|vl/,(x)p)]| (3.7) 
where 

l/T 

*a(x) = ^ f dr^aix^r) (3.8) 

is an imaginary time independent classical field governed by the sombrero 
shaped effective potential V{z) {z = I^E'aP) given by 



V{z) = zr{z) + T V / ^ HDk'' + \ujn\ + r{z)] - ^. (3.9) 



The function r = i/i is to be determined by solving the saddle point constraint 
equation dVjdr = 0, 

Ef dk 1 1 , , 

/ — = -. 3.10 

^^^J 27rDk^ + \cOn\+riz) g 



The scaling limit of equations (13. 9p and (I3.10p can be reached leading to a 
universal, cutoff-independent expression for V{z). First consider Eq. (13.101) 
and note that the T = quantum critical point is at g = gc, where gc is 
determined in the large- A^ limit by 

du; rdk 1 1 ^3_^^^ 



271 J 271 Dk^+ \uj\ gc 
and an ultra-violet (UV) cutoff must be included for finiteness. Defining 



6=^D{ , (3.12) 



cLS cL 

at 




a renormalized tuning parameter, the quantum critical point now resides 
6 = 0, T = 0. Subtracting Eq. (13.111) from (13.101) and performing the sum 
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Fig. 4. The numerical solution of the transcendental saddle point equation (j3.14p 
which will be used in all computations of the effective classical potential ViyZ). The 
symbols were calculated using the approximate solution to r jT found in the metahic 
(Eq. (|3.17|) ) and superconducting (Eq. (|3.21|) ) limits. 

and integral up to a UV frequency cutoff A^^ 



5 + z^D 



dk 

2^ 



In 



. DP 



l^K^ + be + riz)\^<_^be + r(z) 



27cT 



2'kT 



(3.13) 



where if) is the digamma function. In this form, the limit A^ ^ cxd can now be 
safely be taken, and after rescaling to a dimensionless momentum 




■"'I 



+ ^ 1 + 



P + r(z)/T^ 
2^ , 



(3.14) 



This is one of the most important results in the scaling limit, and determines 



r{z)/T as a universal function of S/yT and z\J D/T. A numerical solution of 
Eq. (I3.14P is shown in Fig. HI and we note that it has a minimum possible 
value of — 27r due to the argument of the polygamma function. The effective 
potential V = V{z,6,T) and renormalized mass r = {z,6,T) are actually 
functions of three variables z = I^I/aP, ^ and T. For the sake of brevity, we will 
usually just explicitly indicate their z = I^E'aP dependence whenever possible. 

Substituting the expression for 1/g in Eq. 03.1Up into Eq. 03.90 . and subtracting 
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Fig. 5. The scaling form of the effective potential calculated by including 10^ terms 
in the frequency summation for (5/vT = —0.1,0.0,0.10 (right to left curve). The 
inset shows the convergence properties of the sum for fixed 5/vT and z = |^aP- 



a constant which is independent of z, we obtain using uon = 2TmT 



V{z) 






2n + r{z)/2nT 
^n + r{z)/27TT 



2^ 



(3.15) 



The structure of the effective classical potential indicates that it can be written 
in the scaling form 



T3/2 




(3.16) 



where ^v is a universal dimensionless function. By truncating the sum in 
Eq. (I3.15P at some large, but finite value, the scaling function $v/(^) can be 



evaluated at fixed 6/VT as seen in Fig. [51 For 6/VT = zy D/T = we find 
r = -0.69728 leading to V{0) = 1.5100. 



3.2 Limiting forms of V (z) 



An analytic form for the effective potential can be determined for the limit- 
ing cases 5/\/T — > oo and 6/\/T -^ — oo corresponding to the metallic and 
superconducting phases respectively. 
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3.2.1 Metallic phase 

In the metallic phase, a large positive value of 6 is reached by investigating 
the r{z) -^ oo behavior of Eq. fl3.14p . In this limit, the asymptotic expansion 
of the gamma function can be used, leading to 

r{z) c^ n^ (6 + zy Dj (3.17) 

with the the result plotted using circular symbols in Fig. HI For r{z) ^ 1, 
the effective potential in Eq. (13.151) can be rewritten in terms of a sum over 
Matsubara frequencies. After adding and subtracting the n = term 

yiz) = — 7^\^nz) + ^^}_^\^^^==-2,/UJn]. (3.18) 



In the low temperature limit, the Matsubara summation can be converted into 
an integral 

V{z) ~ -^^) + ^-^[r{z)fl^ (3.19) 

and using Eq. (13.171) we find the temperature independent result 

V{z) ~ iS + z\/D) . (3.20) 



3.2.2 Superconducting phase 

For the strongly ordered phase, 5/y/T -^ — oo and the analysis is less straight- 
forward. As noted previously, r{z) is bounded from below by the first negative 
Matsubara frequency —ui = —27iT, and near this value, the most divergent 
term in Eq. (13.141) produces 

r{z)--ui + T^i6 + z\/D] , (3.21) 

plotted as open squares in Fig. HI With r{z) taking this extreme value, all 
terms in Eq. (I3.15P are well behaved, except the n = 1 term. Extracting the 
culprit 



T 

V{z) = 



\lr{z) + CUi / n=l V \Jr{z) + LOn+l 

(3.22) 
where the sum has been shifted by one. Substituting Eq. (I3.2ip in the first 
term to investigate the divergence, and setting r{z) = —ui in the second, the 
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low temperature form of the interaction potential is given by 

V{z) ~ -27r (5 + z\[d] T. (3.23) 



It seems somewhat surprising that in the superconducting phase, the effective 
potential vanishes linearly with temperature, as T ^ 0. However, at this point, 
we will endeavor to calculate transport only in the quantum critical regime 
where 5/yT <C 1 and can thus safely ignore the irregularity. We will return 
to the issue in Section [33] by investigating the low temperature ordered phase 
through a calculation of the Coleman- Weinberg effective potential at T = 0, 
as well as constructing an effective Ginzburg-Landau potential, near T^. 



3.3 Renormalized classical conductivity 



Now the fruits of our labor become apparent as the dc electrical conductivity 
can be calculated for the effective classical action described by Eq. (13.71) via 



the Kubo formula (1731 ) by reintroducing a real time dependence to the classical 



order parameter and approximating its low frequency dynamics by a Langevin 



equation (15 ll ). The corresponding equation of motion is 



^ = b^2^ - v'm:^,t)mx,t) + n.,() (3,24) 

where / is a complex Gaussian correlated random noise obeying 

(/(x, t)r{x\ t')) = 2T6{x - x')6{t - t'). (3.25) 

By taking a derivative of Eq. (13. 9p and using the saddle point equation (I3.10J1 . 
V'{z) = r{z). Eq. (13. 241) represents the simple Model A dynamics of Ref. (,74 ) 
and should capture the correct quantum critical dynamics whenever the renor- 
malized mass r takes a value such that the Un ^ modes are sufficiently 
gapped. 

The electrical current is defined to be 

J = ie*D (**a^^ - d^^*^) (3.26) 

and thus the uniform dc conductivity can be found from an integral over all 
space and time of the current-current correlation function over the partition 
function Zc 

L 00 

a = - I dx I dt{J{x, t) J(0, 0)) (3.27) 
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Dimensional analysis of the equation of motion above in conjunction with 
Eq. fl3.16p . implies that the classical conductivity obeys the scaling form 



e*2 






where we have inserted the dimensionally correct powers of h and k-^ in the 
final result. This is simply Eq. ( 12. Sp with the replacement given in Eq. ( 12.12^ 
written in terms of our new measure of the distance from criticality 5. 

The scaling function $o-(x) is a smooth function of x through x = and it 
can be determined by finding a numerical solution to the classical equation of 
motion (13.241) for a one-component classical complex field '^{x,t). By employ- 
ing both classical Monte Carlo simulations and a stochastic partial differential 
equation solver, its full time-evolution can be determined. 



We begin by fixing the value of 5/vT and discretize the Hamiltonian described 
by the classical partition function Eq. (13.71) to a unit spaced lattice of L sites. 
In order to solve the full equation of motion, the set of initial conditions, 
given by the equilibrium configurations of the order parameter field \l/ must 
be determined. These can be obtained by classical Monte Carlo methods using 
a large number of different seed configurations where the equilibrium order 
parameters are stored only after a suitable number of Monte Carlo time steps 
(large enough to eliminate any possible autocorrelations) have been performed. 
The set of configurations are then used as the initial [t = 0) states of the 
stochastic equation of motion, Eq. (I3.24p . At each time step, the noise function 
f{x,t) is drawn from a Gaussian distribution and by using the second order 



stochastic Runge Kutta (or Heun) algorithm (J75l ) the time dependence of 
\E'(x, t) is determined. The current-current correlator in Eq. (I3.27P is computed 
as an average over all temporal trajectories of \l/ and the dc conductivity can 
be found after integrating over all space and time. The results are necessarily 
dependent on the size of the spatial and temporal discretization and the final 
value for the scaling function $„- must be finite size scaled in both space and 
time. The resulting value of ^^ directly above the quantum critical point, 
{5 = OT > 0) was found to be $„(0) = 0.07801 ± 0.01. This fully universal 
number is independent of any of the specific details of the particular quasi-one 
dimensional system under consideration. We have also computed the value of 
$cr for a range of 6 near the critical coupling as seen in Fig. El The semi-classical 
result for the physical one-component complex order parameter determined 
here should quite accurately reproduce the real electrical transport in the 
quantum critical regime, and we will use it to benchmark our N = oo results 
in Section [S31 
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Fig. 6. The dc conductivity scaling function $o- near the quantum critical point 
calculated by brute force integrating the current-current correlator measured using 
the Langevin dynamics formalism for a one-component complex field (iV = 1). 
The inset shows the spatial and temporal finite size scaling for a single data point 
corresponding to 6/VT = 0. Note the non-hnear relationship between the spatial 
and temporal mesh sizes. 

3.4 The ordered phase 



To address the physics of the ordered phase we again consider the action 
for an N component complex field \l/a(x,r) with magnitude |\l/a(x,r)| = 1 
but instead of treating the fluctuations semi-classically, we take a a-model 
approach (J76l ). After enforcing the fixed magnitude constraint on \l/a with a 
Lagrange multipher field fi{x, r), its action in the presence of a finite conjugate 
field ha reads 



S 



dx I dT< '^*a{x,T) 



-Ddl + \dr\+i^i{x,T) '^aix,T) - in{x,T 



9[KKi^,r) + K'^ai^,r)] . 



(3.29) 



where we have integrated the kinetic term by parts, and used the abuse of 
notation \dr\ to infer the dissipative \un\ term in frequency space. To derive 
saddle point equations in the large-A^ limit, we explicitly break the 0{N) 
symmetry by choosing only one component of the conjugate field ha to be non- 
zero and equal to h in the N^^ direction. Integrating over A^ — 1 components of 
\l/ and keeping only one component a (not to be confused with the conductivity 
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and dropping the space and imaginary time dependence of fields) 

S,s = {N-1) fdx fdr -a* [-Ddl + \dr\)(y + V (|(t|2, i^) - {ha* + h*a) 

(3.30") 
wliere we liave rescaled the couphng g hy a, factor of A^ — 1. In the hmit 
of large A^, the effective potential V can be determined by the saddle point 
approximation where r = in, 



q ~^J Zn 



along with the constraint equations for r and a 

ar = gh 



(3.31) 



dk 



1 



277 DP + luJ +r 



(3.32) 
(3.33) 



We note that an important distinction between Eq. (I3.32p and (I3.33P and the 
saddle point equation (I3.10p derived in the last section is that here we integrate 
over all Matsubara frequencies. In the absence of an external magnetic field at 
T = 0, the quantum critical point corresponds to the solution a = 0, r = 0, and 
as before defines a critical coupling strength g^ as given in Eq. (13.111) . Again 
we will measure deviations from quantum criticality using the parameter 6 
(Eq. dan). 

Using this definition, the solution to Eq. (I3.32p in zero conjugate field h is 
given by r = and thus from Eq. (13.331) with |(Tp 

9 



\(^o\ 



|o^o 



1-^ 

9c 



6 



(3.34) 



which is clearly only valid in the ordered phase characterized hj 6 < 0. Using 
Eq. (I3.34p . Eq. (I3.33P can be rewritten as 



o" 



\(^o\ 



dk 
'I 2^ 



^Et7 



1 



duj 



(3.35) 



^.^^ Dk^ + \ujn\ + r J 2tt Dk^ + \ujr, 
and by a method identical to the one used when integrating over all Un ^ 0: 



O" 



9 



^5+gf 
n J 



dk 

2^ 



Dk"^ + r' 



V' I 1 + ~ : ^ • I +\n(^\ 



-kT 



27iT 



\DkV Dk"^ 



(3.36) 



where ip is the digamma function. This expression can be inverted numerically 
to provide r as a function of |crp and 6. After this has been accomplished, 
Eqs. (13.311) and (13.331) can be combined to give the finite temperature effective 
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potential 



.,|.p,.,^rE/|[,n(i.-: 



r{\a\\6) 



r{\a\\5) 



(3.37) 
where a constant term independent of r has been discarded. 



3.4-1 Zero temperature effective potential 



At zero temperature, the frequency and momentum integrals in Eq. fl3.35p can 
be written in an isotropic fashion. For a finite conjugate field, both |crp and r 
are nonzero, and dropping the explicit |crp and 5 dependence of r 



V{\a\\5) 



d^p 



^3/2 



In 1 + - - -^—- 
\ p"^ I p'^ + r 



SnVD 
where the saddle point equation fl3.33p can now be solved as 

9 



(3.38) 



0" 



koP + 



TiVD 



(3.39) 



which indicates that a solution exists only for |o"p > |crop, requiring the pres- 
ence of a non-zero conjugate field h. In the ordered phase, this equation cannot 
be solved for |o"p < |crop, and hence the effective potential defined below will 
not be valid in the weakly ordered regime. 

The zero temperature Coleman- Weinberg effective action for a quantum field 
^ is defined to be (1771 ) 



r[^e/] = -S,^[^>a]- j dx j dT{h*m^i + hmli 



(3.40) 



such that it is the function whose minimum gives exactly ^^ = (^)- To lowest 
order in perturbation theory it is simply the classical potential energy, but is 
modified by quantum corrections at higher order. Using Eqs. (I3.30p . (13.380 
and (]3.39p . we obtain 



K 



e// 



fi(iv-i) sVd 



cr 



(3.41) 



where Vt denotes the system volume in space-time. As noted above, this effec- 
tive potential is only valid for |crp > |croP, as it has a minima at |crp = (for 
(5 < 0) and a point of inflection at |(jp = |o"op. 
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Fig. 7. The shaded portion shows the region of integration after the implementation 
of a symmetric infrared cutoff in frequency and momentum. 

In order to find a solution for |crp < |o"op, and to derive a Ginzburg-Landau ef- 
fective potential for the description of any slow degrees of freedom, we should 
not integrate over all degrees of freedom, but only over those with a wave- 
length smaller than some cutoff A~^. This is necessary due to the fact that 
when integrating over all Matsubara frequencies we are restricted to r > 



and thus can never access |fTp < |o"oP- In Ref. ( 1l2l ) a similar viewpoint was 



expressed, and the cutoff was taken to be of the order of the zero temperature 
superconducting coherence length. We will follow this approach here with the 
addition of a step that ensures self-consistency, and the cutoff is implemented 
symmetrically in momentum and frequency space, consistent with dynami- 
cally critical exponent z = 2. A schematic diagram showing the shaded region 
of integration can be seen in Fig. [71 The zero temperature effective potential, 
with r < is now defined as as 



K 



eff 



D J Wr 

\/d\ 



\ p^ j p2 



9 



(3.42) 



Differentiation with respect to r gives rise to the modified saddle point equa- 
tion (for negative r) 




(3.43) 



where |o"op is now modified from Eq. ( 13.34^ as a result of the infrared cutoff 

koP = -4i + ^- (3.44) 



/d tt 



2 ■ 



In order to gain intuition about the relative size of the cutoff and the effective 
mass (Lagrange multiplier) r, we calculate the zero temperature supercon- 
ducting coherence length ,^(0) as a function of r. It is determined as usual by 
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the relation 



dV, 



eff 



d\a\ 



k|2=0 



D 1 

Taoy 



(3.45) 



By explicitly differentiating the effective potential Eq. (13.421) and using Eq. fl3.43p . 
we find 



dV, 



eff 



d\a\ 



k|2=0 



1 

-i 
9 



-ri a 



0), 



and combining the last two equations leads to the relation 



(3.46) 



m 



D 



\|r(|aP = 0)r 



(3.47) 



Using this definition at T = and | a p = we can now determine the coherence 
length ,^(0) self-consistently from Eq. (I3.43P 



m 



TT^\6\ 



In 






2Am 



(3.48) 



where Eq. (13.441) has been used. Note that this equation has a solution for all 
choices of A such that 1 < ^(0)A < 6/5 with ,^(0) —>■ cxd logarithmically as 
A^(0) -^ 1 and ^(0) ^ as A^(0) -^ 6/5. Parameterizing A^(0) = 1 + e where 
e <^ 1, and defining 



m = -. 



TT" 



In ( 1 + - 



2(1 + el 



(3.49) 



the zero temperature coherence length is 



m 



1^1 



/w- 



(3.50) 



Due to the logarithmic divergence as e — > 0, one possible choice of e ~ 
1.4 X 10~^ gives /(e) ~ 1 leading to the simple relation ^(0) = \ D/\5\ or 
A c:^ \5\/vD. The preceding arguments now allow us to express the effective 
potential (Eq. (13.421) ) in terms of |cr(r < 0)^ < |crop and 5. If we choose e such 
that /(e) = TT^, i.e. ^(0) = v^7rVl<^| then the saddle point equation ([333]) 
simplifies to 



cr 



1 + 



vr^ 



1^1 



w |r|arctanh 




(3.51) 



which when solved for numerically for r can be substituted into Eq. (I3.42p to 
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Fig. 8. The effective quantum potential at T = calculated via the Coleman- Wein- 
berg procedure for \a\^ > |croP ^-nd through the self-consistent insertion of an in- 
frared cutoff for |crp < |o"o|^ with /(e) = tt^ for 6 = 0.0, —0.40,0.60 (inner to outer 
curve) . 

give the effective potential 



K 



eff 



+ 



37r2 



cr 



9 



TT" 



-r5 + (5^ In f 1 + ^ ) + Vr^\n 




(3.52) 



This result, valid at T = and |o"p < |o"oP can now be combined with 
Eq. (13.411) which is valid for |o"p > |o"oP to obtain the effective potential ev- 
erywhere at T = and 6 < corresponding to the ordered or superconducting 
state, as seen in Fig. [HI 



3-4-2 Finite temperature Ginzburg- Landau potential 

Having computed the form of the infrared momentum cutoff A self-consistently 
as a function of 5 and after investigating the form of the effective potential at 
zero temperature, we now move to finite temperatures and consider expanding 
around some critical temperature Tc for the ordered phase. The usual form for 
the potential is posited 



VGL = Vo + ao{T-T,)\a\ 



Ipw\^ 



+ ■ 



(3.53) 
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Fig. 9. The Ginzburg-Landau potential of Eq. (|3.53|) for 5 = 0.0, 
to top curve) using the values of Tc, oq and (3 found for /(e) = vr 



0.4,-0.6 (bottom 
2 at T = O.BT^. 



and we will endeavor to evaluate T^., ao and (3 in terms of the parameters g, D 
and 5. The goal of such a procedure is to derive an effective Ginzburg-Landau 
theory for the superconducting state near T^ with quantum renormalized co- 
efficients. By multiplying this potential by the finite temperature Ginzburg- 
Landau coherence length ^(T), an effective free energy is found from which the 
barrier height for a thermally activated phase slip can be determined directly 
from the LAMH theory. The details of the calculation of the renormalized cou- 
pling constants is rather involved and has been included in Appendix O The 
resulting dimensionless interaction potential is shown in Fig. [9]for T = CST^. 
We find a relatively steep double well potential which unsurprisingly has both 
a barrier height and order parameter expectation value that depends on the 
distance from the critical point 5. 



3.4-3 Quantum renormalized LAMH theory 

Having derived a effective Ginzburg-Landau potential (with engineering di- 
mensions of energy divided by length) in Appendix [Cj we may convert it into 
a free energy functional by multiplying by the finite temperature coherence 
length 

e(T) = m (l - ^) • (3.54) 

Due to the presence of phase slips below T^, phase coherence of the wire is lost 
and it necessarily breaks up into L/^{T) independent segments, that interact 
via Josephson coupling near the finite temperature transition. After discarding 
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a constant, the functional is given by 

^^^^^> = "° ^^ I' - tJ '^' ~ ^^?^ I' " tJ i^i 

(3.55) 
with a singular temperature dependence in the quartic term. We restrict the 
analysis to temperatures that place the system in the LAMH region of the 
phase diagram displayed in Fig. [3l such that 1 — T/Tc ^ 1. The rescaled 
distance from the critical point is negative {6 < 0), and the dimensionless 
coefficients are 



ao = gao (3.56) 

^(3 (3.57) 

where ao and /3 are given in Eqs. (10.50 and (IC.lip respectively. 

The free energy barrier height for a single phase slip event can now be evalu- 
ated using the coefficients computed directly from the non-linear sigma model 
version of the full quantum theory, 

2 7i72 f(AT3/-2 / T\3/2 



a a 



■lf{^)T!"(, T\ 



AF = -^ = -^ ^^^/ 1 (3.58) 

which can be written in terms of a dimensionless scaling function of two unique 
scaling variables, the first expressing the classical and the second the quantum 
nature of the renormalized barrier height 



AF = T^^^F h^, ^^ • (3.59) 




AF can now be directly inserted into the LAMH theory in place of the phase 



slip barrier height calculated by Langer and Ambegaokar (11 ll ). with the rest 



of the arguments leading to the LAMH resistance remaining unaffected. 

Before doing so, we briefly remark on the seemingly surreptitious form of 
Eq. ( 13.59p . The ability to write the barrier height as a cutoff independent 
scaling function with no explicit dependence on A or /(e) is a consequence of 
the fact that by flxing the dimensionless variable 5/ yfTc a unique value of e and 
thus /(e) can be found from Eq (1C.3I) . This is demonstrated in Fig. [TOl where 
it appears that $af is nearly a linear function of (5/v^, with a slight kink 
for T = 0.8Tc. Therefore, using the same method as discussed in Section [TTT] 
we can write the quantum renormalized LAMH conductivity in terms of the 
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Fig. 10. The dimensionless free energy barrier height corresponding to Eq. (j3.55p as 
a function of a single scahng variable 5/yjT^. 

scaling function $af as 



e2 , 1 

O-QRLAMH = 47r— ^0) 



X exp 



$AF I — 




j^s 3/2 



T 



-3/2 



(3.60) 



or 



ctqrlamh 



-m $ 



QRLAMH 




(3.61) 



where $qrlamh is shown in Fig. [H] for three values of bj \fTc- The arguments 
discussed in Section IB. II provide a recipe to convert the parameter b to the 
physical pairbreaking frequency a for a given experimental geometry. The sup- 
pression of the critical temperature as a result of an external magnetic field 
directed parallel to the wire has already been observed in Ref. (1171 ). Thus, 
in principle, the relationship between 8 and Tc could be determined experi- 
mentally from a fitting procedure, and the expression for the conductivity in 
Eq. (I3.6ip could be applied to the experimental transport results with one less 
fitting parameter than the form currently used in Eq. (II. 3p . 

A striking difference between the usual LAMH form of the conductivity and 
the quantum renormalized expression in Eq. (I3.6ip is the absence of any ex- 
plicit dependence the number of transverse channels, N ^ in the latter. This can 
be most easily seen by using Eqs. (I3.50p and (10.40 to write the dimensionless 
variable (5/v^ in terms of Z^/^^(0)Tc, the natural energy ratio suitable to our 
analysis. Due to their dependence on different microscopic length scales, it is 
somewhat difficult to quantitatively compare ctlamh with its quantum renor- 
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Fig. 11. The quantum renormalized LAMH scaling function for the conduc- 
tivity plotted as a function of the reduced temperature for three values of 
6/^/T^ = -1.25, -0.50, -0.25 (top curve to bottom). 

malized value ctqrlamh- Qualitatively, the inclusion of quantum fluctuations 
leads to a softening of the free energy barrier resulting in an enhanced phase 
slip rate consistent with recent calculations made by Golubev and Zaikin (1781). 



3.5 Large-N expansion 



Moving away from the ordered phase, at low temperatures, at a distance from 
the quantum critical point {VT <^ 5), quantum fluctuations are large and 
there will be a flnite number of Matsubara frequencies which lie below the 
crossover energy. Now, the classical model of Section 13.11 is no longer adequate 
to describe the contribution of superconducting fluctuations to transport. In 
this regime of the phase diagram however, a direct 1/N expansion on the full 
quantum theory can be attempted. Starting with Eq. (13.11) and decoupling 
the quartic interaction with a Hubbard- Stratonovich fleld /i we arrive at the 
effective action 



dx / dr 



D\d^'^a{x,T)\^ +ifl{x,T)\'$a{x,r)\^ + —fi\x,T) 

2u 



.a 



+ i—fi(x,T) 
u 



Tj2jdx\uJn\\^a{ 



(3.62) 



Integrating out \l/a over its now quadratic action in the partition function 
Z = Tt exp(— iS), we as usual recognize an overall factor of A^ which allows 
the functional integral over the Hubbard-Stratonovich fleld /x to be performed 
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in the saddle point approximation where we identify R = ifi. In the universal 
limit, the new quadratic effective action is given by 



SR = TY^J—{De + \u:^\+R)\^,{k,u;n)\' 



where the 'mass' R is defined by the saddle point condition 

1 rdk_.^ 1 

9 



—TY^ 



(3.63) 



(3.64) 



Unl + R 

The evaluation of R is straightforward and follows our derivation of Eq (13.141) 
6 r dk 



T 



27r2 



, , 27r\ , / k'^ + R/r 



TX 



A;2 + R/T 



(3.65) 



In general, this expression must be inverted numerically to determine R/T as 
a function of 5/\/T as is shown in Fig. [T21 However, at the quantum critical 
point (QC, 5 = 0) and in the metallic (M,(5 -^ cxd) and superconducting (SC, 
5 —>■ — oo) we can analyze Eq. (13.651) along the same lines as was done for 
Eq. fl37[4D in Section O We write R = T^r{6/VT) and find the following 
results 

1/4x2 ; 

(3.66) 



^r{x) 



X — > — oo 
0.625 ; X < 1 



2 2 

TV X ; X 



CXD 



Understanding the behavior of the effective mass will be crucial in the follow- 
ing sections where we first define and then calculate the thermal and electrical 
transport coefficients in the zero frequency limit. 



3.5.1 Thermoelectric transport 

The electrical (cr) and thermal (k) conductivities and the Peltier coefficient 
(a) are defined in terms of the electrical Jq and thermal J2 current densities 
via the relation 

^j,j [aTk)[-VT 

where E is an external electric field, VT is an imposed temperature gradient 
and k/T = k/T — o? jo. The current operators can also be defined in terms 
of the derivatives of our large- iV action iS/j 



(3.67) 



Jo 

J2 



Mo 
QSk 



(3.68) 
(3.69) 
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Fig. 12. The renormalized mass R plotted as a function of the rescaled distance 
from criticality 5/^/T. The symbols refer to the analytic results of Eq. (j3.66p in the 
metallic and superconducting limits. 

where Aq is the scalar electric potential and A2 is the thermal vector potential, 
after we have made it gauge- covariant through the introduction of these two 
gauge fields (J79l ) via the replacement 



leading to 



dx ^V = d^ ~ie*Ao{x,T) - iA2{x,T){idr] 



Jo = ^e*D [rVi^ - ^(^^)*] 
j2 = D [dripiV^py + Vijd^ij*] 



(3.70) 



(3.71) 
(3.72) 



where e* = 2e is the charge of a Cooper pair and we have specialized to 1 + 1 
dimensions. 



The quantum Kubo formula (l73l : ISOl ) can be used to obtain results for the 
thermoelectric conductivities at external complex frequency iuon (where we 
ignore the Peltier coefficient as its dc part will turn out to be identically zero) 

m 



QpiiUn) 



d I dS 



LUnTP dAp \dAp^ 



Ao=A2=0 



1 

'uj.TP 



p 



(3.73) 
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where the currents are defined by 



(3.74) 
and p = corresponds to the electrical conductivity while p = 2 defines the 
thermal conductivity, i.e. a^iun) = Go{iuJn) and K{iujn)/T = Q2{i^n)- 



The conductivities are more easily expressed in terms of a one-loop polariza- 
tion function, 



4L)2e*2-p 

Qp{iuJn) = ——ICp{iuJn) 

which contains both paramagnetic and diamagnetic contributions 



(3.75) 



ICpiiUn. 





fen 



1 
(|cu„|+£)P + i?)2j 



{\ujn\ + Dk^ + R){\uJn + en\+Dk^ + R) 

(3.76) 



where a solid line represents the bare propagator Go(fc,co'n) = {Dk"^ + |a;„| + 
R)~^-i an open circle corresponds to a term linear in the potential Ap and an 
open square to a term quadratic in Ap. Employing the spectral representation 
of the bare propagator 



A{k^uj) = —2lm.Go{k, iun -^ u + irj) 

2uj 



cj2 + (Dk^ + Rf 



(3.77) 



where a ItUn] dependence along the imaginary frequency axis becomes —iuj just 
above the real frequency axis. The polarization function is then given by 



.,,.,^ri:/|.vt/ 



dk 2 [ duJi r duj2 fuJi + uJ2\^ 
27r J ^ J ^ V 2" 
1 



A{k,uJi)A{k,uj2) 
1 



(ie„, - iui) [i{en + iUn) - ^^2] («e„ - u;i) (ie„ - UJ2) J 

(3.7; 



where we have made the replacement (e„ -|- u;„/2) -^ (tui + co'2)/2, due to the 
temporal non- locality of J2 (l79l ; ISlf l. Performing the Matsubara summation. 
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and analytically continuing to real frequencies yields 



'=^(-'^)-/S/^/^(^r-^(^-)-^(^-'-)^-= 



27r V 
X [n{uji) -n{uJ2)] I 



1 



1 



(3.79) 



UJ2 — UJi — UJ — trj UJ2 — UJiJ 

where n(u;) = (e^^^ — 1)^^ is the Bose distribution function, and r/ is a positive 
infinitesimal. After checking that the delta-function contribution to Re^p(cc;) 
at zero frequency is proportional to the external frequency, (it vanishes as 
uj —^ 0), we can combine Eqs. (I3.75P and 03.791) to give the remaining regular 
part 



ReQp(Lu) 



4^2g*2-p ^^ ^^^^^ _ ^^^ ^ ^^^ 



rpp 



TT 



UO 



n 



2, 



X 



dk 



k^n{n + uj) 



271 [fi2 + p^2 + i?)2][(^ + ^)2 + p^2 + f>y 



(3.80) 



T/o;, 



The classical limit of Eq. (I3.80p can be found by replacing n{uj) 

but here we directly perform the limit uj ^ and obtain the quantum dc 

conductivities 



Re^„ 



D'e 



2.*2-P 



dQ 



n^+p 
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8TP^ 



271 



sinh {n/2T) ^^2 + ^2 (^^ + VTFT^R^ 



3/2' 



(3.81) 



which is the major result of this section. From this expression it is immediately 
clear that the Peltier coefficient a (corresponding to p = 1) is identically zero 
by symmetry, and thus k = k. 

In addition to the scaling function $o- defined for the electrical conductivity 
in Eq. 03.28P and given by 



-,*2 



hD 



:$. 



6 



the thermal conductivity must obey a similar form 



(3.82) 



T 



A-2 



hD 



h \ k^T 



$. 



5 



^/hk-QT ) ' 
where we have re-inserted the appropriate factors of h and k-Q for clarity. 



(3.83) 
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Fig. 13. The solid lines show the N = oo universal scaling functions the electrical 
(p = 0, top curve, left axis) and thermal (p = 2, bottom curve, right axis) conductiv- 
ity calculated by integration of Eq. ()3.8ip . The symbols show the effective classical 
scaling function for the electrical conductivity calculated in the Langevin formalism 
in Section 13.31 and previously shown in Fig. [6] for a one component complex field. 



The bjyT dependence of ^^ and $k can be found by numerically inverting 
Eq. (I3.65p . (Fig. [T2|) and the result is shown in Fig. [131 Both are smooth func- 
tions of 5 /yT throughout the quantum critical regime with the dimensionless 
value of the electrical conductivity being two orders of magnitude larger than 
the thermal conductivity right above the quantum critical point. The large- A^ 
electrical transport can also be compared with the previously calculated value 
found in Section [331 for A^ = 1. There is good agreement between the two, and 
both results have the same 5 /yT dependence near 5 = indicating that the 
correct physics are manifest even aX N = oo. With the full numerical scaling 
functions computed, the temperature dependence of the thermoelectric trans- 
port can be determined by fixing 5 leading to the results displayed in Fig. [H] 
and [151 The singular correction to the electrical conductivity clearly shows be- 
havior consistent with a quantum phase transition between a superconductor 
(diverging conductivity) and metal (finite or vanishing conductivity). More- 
over, for a fixed but small positive value of 5, placing the system in the quan- 
tum critical region of Fig. ([2]), (inset Fig. [H]) as the temperature is lowered, 
a clear signature of crossover behavior can be seen from the non-monotonic 
temperature dependence of the conductivity. As the system leaves the strongly 
fluctuating quantum critical regime for the low temperature metallic regime, 
AL contributions from Cooper pairs are suppressed and their positive con- 
tribution to the conductivity disappears. As mentioned in the introduction, 
experiments (1821 : l33l ) have already seen evidence of such non-monotonic re- 
sistance in the metallic regime, strongly supporting the crossover picture we 
present here. The exact form of the temperature-dependent crossovers can 
be determined by further investigating the limiting forms of the two scaling 
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Fig. 14. The temperature dependence of the rescaled dc electrical conductivity at 
fixed values of <5 = -0.50,-0.20,-0.10,0.00,0.10,0.20,0.50 increasing from top to 
bottom. The inset shows behavior non-monotonic in temperature near criticality. 




Fig. 15. The temperature dependence of the rescaled dc thermal conductivity at 
fixed values of 5 = -0.50,-0.20,-0.10,0.00,0.10,0.20,0.50 increasing from top to 
bottom. The inset shows behavior non-monotonic in temperature near criticality. 

functions $o- and $k- 



3.5.2 Finite temperature crossovers 

Analytical forms for the electrical and thermal conductivity can be found in 
three distinct limits using the results of Eq. (13.661) . The first is deep in the 
superconducting (SC) regime where S/y/T -^ — oo or R/T -^ 0. The next is 
at the quantum critical point (QC), where 6 = or R/T = 0.624798 and the 
final regime is on the metallic side of the transition (M) where S, R/T -^ oo. 
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SC 


QC 


M 


a 

KJT 


i/r2 


1/Vt 
1/Vt 


T2 
T2 



Table 2 

A summary of the temperature dependence of the electrical and thermal conductiv- 
ity in the superconducting (SC), quantum critical (QC) and metallic (M) regimes. 

For low temperatures, Eq. (I3.81|) can be evaluated in these three limits leading 
to the approximate analytic scaling behavior 



$^(x) 



^Ax\ 



x' 


; SC 


0.217997- ■■ 


■QC 


(127r4)"ix-5 


; M 


4v^^ [2) 


; SC 


0.24592- ■• 


■QC 


(157r2)-ix-5 


;M 



(3.84) 



(3.85) 



which can be used in Eqs. (I3.82p and (13.831) with x = d/^/hk-oT. The leading 
order temperature dependence of these results is summarized in Table [H The 
finite temperature crossover behavior near the z = 2 SMT is characterized by 
the conductivity increasing like 1/T^/^ at high temperatures while the system 
is in the quantum critical regime of Fig. [3] and finally decreasing as T^ after the 
system has fully returned to metallic behavior. Although we are about to show 
that a microscopic theory reproduces the T^ metallic conduction, the 1/vT 
dependence of the conductivity is not present in the simple Gaussian theory 
and an accurate determination of the full crossover phase diagram necessitates 
the inclusion of interactions between Cooper pairs. 

We now comment on a shared regime of validity between our large- A^ theory 
and the disordered electron perturbation theory of Ref. (l43l ). A closer inves- 
tigation of the dc electrical conductivity in the metallic regime with careful 
attention to all prefactor yields 



a 



T 12/?5/2 • 



(3.86) 



We note that upon comparing the large-A^ propagator of Eq. (I3.64p with 
Eq. (4) of Ref. (1431 ) that the mass R used here is exactly double the mass 
a employed by Lopatin et al, i.e. R = 2[a — ac{T)]. Having made this identifi- 
cation, we may compare Eq. (I3.86P above, with the finite temperature fluctu- 
ation correction to the normal state conductivity computed via diagrammatic 
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perturbation theory (Eq. (8) in Ref. (J43|)), and find exact agreement. After 



a rattier lengtliy calculation it was confirmed that perfect correspondence is 



also found for the thermal conductivity in this limit (16 ll ). The concurrence 
between the two theoretical approaches in this limit is a result of an approxi- 
mation made in the diagrammatic calculation involving an infinite sum over a 
class of ladder diagrams which turns out to be equivalent to the large- A^ limit 
taken here. 



3. 6 Wiedemann- Franz ratio 



Having computed the temperature dependence of both the thermal and electri- 
cal dc conductivity via a large- A^ expansion, it is an opportune time to consider 
what these results may indicate about the physics of the SMT. This is mostly 
easily accomplished through an analysis of the validity of the Wiedemann- 
Franz (WF) law, which states that the low temperature limit of the ratio 

W^^ (3.87) 

of the thermal and electrical conductivities of metals if finite and given by the 
universal Lorenz number 

,-^^)\ (3.88, 

This remarkable value relates macroscopic transport properties which can be 
found from an analysis of current-current correlation functions to fundamental 
constants of nature. The prefactor vr^/S is fixed solely by the Fermi statistics 
and charge of the elementary quasiparticle excitations of the metal. Physi- 
cally, any constant value of W indicates that elastic collisions dominate as 
all scattering events are necessarily charge conserving. Both the temperature 
independence of W as well as the actual value of Iq has been experimentally 
verified to high precision in a wide range of metals (83-), and realizes a sensitive 
macroscopic test of the quantum statistics of the charge carriers. 

In addition to simple metals, the value of the Wiedemann- Franz ratio is known 
in some other important strongly interacting quantum systems. In supercon- 
ductors, which possess low energy bosonic quasiparticle excitations, a is in- 
finite for a range of T > 0, while k is finite in the presence of impurities 



(1841). and thus IV = 0. At quantum phase transitions described by relativistic 
field theories, such as the superfluid-insulator transition in the Bose Hubbard 
model, the low energy excitations are strongly coupled and quasiparticles are 
not well defined; in such theories the conservation of the relativistic stress- 



energy tensor implies that k is infinite, and soW = oo (l85l). In other words, 
any quantum critical point which exhibits Lorentz or Galilean invariance will 
have an infinite thermal conductivity since the boosted thermal distribution 
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will never decay (15 ll ). 



Catelani and Aleiner (l86l : 1871 ) have recently investigated interaction corrections 
to the Lorenz number in disordered metals. They found that neutral bosonic 
soft modes resulting from interacting electron-hole pairs contribute to the 
transport of energy leading to temperature dependent deviations from Iq, with 



the corrections being larger in lower dimensions. Li and Orignac (l88l ) computed 
ly in a disordered Luttinger liquid, and found deviations from Iq leading to a 
non-zero universal value at the metal-insulator transition for spinless fermions. 
Finally, Fazio et al. have computed the effects of plasmon scattering on the 



Lorenz number of thin wires coupled to reservoirs (J89l ) where charge-energy 



separation leads to a violation of the Wiedemann- Franz law. 

For the SMT considered here, an examination of Table [2] immediately leads to 
the observation that the Wiedemann-Franz ratio is temperature independent 
in both the quantum critical and metallic regimes. Remarkably, all important 
couplings between bosons and fermions scale to universal values, and conse- 
quently, by studying the analytic form of a and k/T given in Eqs. fl3.82p and 
we find the universal constant 




PFqc = (0.28203- ••) (^J (3.89) 

in the quantum critical regime whereas in the metallic region of the phase 
diagram 

(3.90) 

Both of these corrections are smaller than the Lorenz number and thus it 
appears that the Cooper pairs tend to carry more charge than heat. 

These results for the Wiedemann-Franz ratio conclude this section, but in the 
next section the theory presented here for N = oois extended to the first order 
in 1/A^. We will exploit the anomalous scaling dimension of the dynamical 




critical exponent z to find additional universal corrections to Eq. (I3.89p . 



4 1/A^ corrections 



This section is quite heavy on calculational details and can be skipped by the 
casual reader. The main results include the derivation of a critical theory for 
a finite A^ component complex field '^a governing the fluctuations of Cooper 
pairs near a superconductor-metal transition. This theory is used to systemati- 
cally compute the 1/A^ corrections to critical exponents and the zero frequency 
transport coefficients calculated in the previous section at A^ = oo when the 
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coupling parameter which drives the SMT attains its critical value. We will 
find (Eq. fl4.79p ) that although the individual values of the electrical and ther- 
mal conductivities are not universal, instead depending depend explicitly on 
an ultra-violet cutoff, their Wiedemann- Franz ratio is a pure, temperature in- 
dependent universal number that characterizes the most singular corrections 
to transport in the strongly fiuctuating quantum critical regime 

4-1 The critical theory 



We begin by reintroducing the strong-coupling effective action of Eq. fl3.4p for 
an A^-component Cooper pair operator \l/a 



S9 = -TJ: / ^{k' + |..„|)|vl/„(fc,a;„)|2 (4.2) 

q t-^J 27r 



where distances have been rescaled by a factor of the square root of the ef- 
fective diffusion constant D and "hard spin" constraint |\l/(a;, r)p = 1 must 
be enforced. Imposing the delta-function constraint via a Lagrange multiplier 
/i and performing a rescaling of the field \l/a — *> ^/g'^a leads to the partition 
function 

X — dx dr 



*:(x.r)(-a^ + |S.|+.,(x.r))*„(..r)-^.M..r) 



(4.3) 



where we have again used the notation \dr\ to infer \ujn\ after Fourier trans- 
forming. Integrating out the \l/a fields, 

Z = I P/i exp < — A^ Tr In i—dl -|- |9t-| + 2/i(x, r) j I dx I dT^{x, r) 

1.4) 

and as done previously, for A^ large, we can approximate the functional integral 
over /i by its saddle point value defined to be r = ifi leading to 

I r dk 1 

^" 1;^] 2^P + |^„|+r" ^^-^^ 

This is an auspicious point to make a brief comment regarding the relation- 
ship between the notation introduced here and that of the previous section. 
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Eq. (14. 5 p is identical to Eq. (13.641) with the replacement oi r ^ R. However, 
in this section, unless otherwise stated, we will be considering a critical the- 
ory (whether at zero or finite temperatures) with the coupling g equal to its 
critical value g^ which will be shifted from its A^ = oo value by a correction 
of order 1/A^. As a result, the effective mass r will be also corrected from its 
N = oo saddle point value. This will be made more explicit soon, but for now 
we simply indicate that r = R + 0{1/N) with R equal to its A^ = cxd value 
defined by Eq. (13.641) with g = gc- 

Let us now focus on fiuctuations around the saddle point by defining ifi = 
r + iX, and after expanding to quadratic order in A and noticing that with the 
help of Eq. (14. 5 p all linear terms cancel we have 



-A^ 



Z = j T>\ exp 
where 






2tt [{k + g)2 + \ujn + en\ + r](g2 + |e„| + r) 



(4.6) 
(4.7) 



can be thought of as the propagator for a A field leading to 1/A^ fiuctuations. 
Upon examination of Eq. (14. 6p . it is apparent that we could have simply 
started from a partition function for the original field "^a with an additional 
interaction term such that its diagrammatic expansion is equivalent to that of 
Eq. (HD, i.e. 



+ iX{x, T)\^aix, t)P H / dx' / dT'X{x^ t)I[t{x — x',t — r', r)A(x', r') > 



leading to the effective action in momentum space 



Sr=TY.j 



dk 

2^ 



{k^ +\uJn\+r)\^a{k,UJn)? + —\X{k,UOn)\^YiT{k,UJn,r) 



dq 



+ ^E / ^^^(^' ^«)^a(?> en)A(fc - g, ^n - Cn) 



(4.9) 



where we note that in order to avoid double-counting, the A or fluctuation 
propagator 11^ cannot have a self-energy contribution of a single \l/a bubble 
(since it has already been included). Thus, performing a direct 1/A^ expansion 
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from Z for G{k, oOn) = {\'^a{k,ujn)\'^) Q we have (to order 1/A^) 



d 



G{k,UJn) 



+ 



(4.10) 



where a sohd hne is equal to (A;^ 



UJr, 



+ r) ^, a dashed hne equal to Ut/N 



and a solid dot represents the interaction vertex i. There is no tadpole graph 
as it is already included in the 1/A^ correction to the effective mass r. The 
third graph has two loops, but is only of order 1/A^ as any closed \l/a loop 
gives a factor of A^. Combining these graphs leads to the full expression 
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27r 
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{p'^ + Wn\+ r)2[(j9 + g)2 + |e^ + z/^l + r] ' 



(4.11) 



4-2 Quantum critical point 



At T = 0, the critical point Qc is determined by the condition G ^(0,0) = 0, 
r = Vc- Keeping terms only up to order 1/A^ 
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dk r duj 



1 1 

+ 
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47rJF + 2|cu| 



2J|cjK/F + 2|cu| -A;2-3|cu| 



(4.13) 



with details given in Appendix [Dl Note that IIo (0,0,0) is infrared divergent, 
but this will shortly cancel out of observable quantities. Inserting the expan- 
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sion for Vc in Eq. (14. 5p . we obtain 

1 r dk r duo 1 

Jc~ J ^ J 271 k^ \uj\ ' N J 2n J 2n 

X 



1 r dk r duj 



dq f de 



1 



2ti J 27rno(g,e,0)(A;2 + |cj|)2 



which is free of infrared divergences. 



4- -3 Quantum critical propagator 



g2 + |e| {k + g)2 + |ti; + e| 

(4.14) 



Moving to finite temperatures, but setting g = g^, we write r = R + Ri with 
-Ri ~ 0{1/N). As mentioned previously, R is determined by setting r = i?, in 
Eq. (14.51) when g = gc takes its A^ = oo value 



^E 



dk 



2nk^+ \un\ + R 



dk r duj 



271 J 27rA;2+ \u\ 



(4.15) 



We have seen this equation before in Eq. (I3.64p with g = gc and can thus 
express it as Eq. (13.651) with 6 = giving an equation that can be inverted to 
uniquely determine R/T, 







dk 

2^ 



-T _,Ji + *^ + «- 



In 



F + i? " V~ ' 27rT I ' '"\27vT . 
where iplx) is the polygamma function. Solving numerically we find 

- ~ 0.624798. 
T 



(4.16) 



(4.17) 



Returning to Eq. (14. 5 p we can write (to order 1/A^) 

1 „^ r dk 1 

9c 



TE 



27rfc2+ \uj^\+R + R^ 
^ J 2nk^+ \uJn\ + R 



(4.18) 



which can be compared with our expression for 1/gc in Eq. (I4.14p order by 
order to yield 
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(4.19) 
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and note that vT 11^(0, 0, R) is a finite universal number given by (see Ap- 
pendix [Dl) 



TUTiO,0,R) 



4(27r)3/2 



Ck 



3 R \ /3 /? 



+ CU 



+ 1 



(4.20) 



.2' 2txT) ' ^ V2'27rT 
where C(^™5 2;) is the Hurwitz Zeta function. Inserting everything in Eq. (14.111) 

G-\k,UJ„) = k'^ + \uJr^\+R + Rl + ^{k,UJn) (4.21) 

where the self energy T,{k,ujn) is defined to be 



Ti{k,UJr, 



AT 2^ 



N^J 2TiIiT{q,en,R)[{k + qf + \uJn + en\+R q^ + \en\+R 

(4.22; 
such that S(0, 0) = 0, and 



Ri 



dq f de 



dk r duo 



A^nT(0,0,i?) [ J 271 J 27rno(g,e,0) y 271 J 27r (P + |cj|)2 
1 1 



+ TE 



g2 + |e| (A; + g)2 + |cj + e| 
dq 1 ^tr^ C dk 



271 UT{q,e,R) 
1 



tj: 



27r(A;2 + |w„|+i?)2 
1 



{q^+\en\+R) {k + qy + \uJn + tnl + R) 



(4.23) 



This is equivalent to Eq. (4.6) in Ref. (J90l ) . Now, both the inner integral or 
sum over {k,uj) and (A;, a;„) is ultraviolet convergent, but the outer integral 
appears to be divergent. By first introducing an ultraviolet momentum cutoff 
and evaluating the integrals numerically using adaptive mesh techniques for 
fixed A the UV limit can be investigated with the result that Ri/T converges 
to the universal finite value 



Ri 0.1069 

T ~ N ' 



(4.24) 



Therefore, as described in Ref. ( l65l ) and Section 12.11 as a consequence of the 
scaling relation z = 2 — 7] where r] is the anomalous dimension of \l/a, the 
uniform static order parameter susceptibility 



X 



dx dT{^:{x,T)^MO)) 



(4.25) 



is determined by the value of fc^T alone. Using Eqs. (14.171) . (I4.2ip and (14.241) 

we find 

0.107^ 



X 



kBT 0.6248 + 



A^ 



(4.26) 
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4-4 Critical exponents 



With a quantum critical theory firmly established, we may now investigate 
any possible 1/A^ corrections to the large- A^ critical behavior characterized by 
exponents z = 2 and u = 1. Such corrections can be obtained by exploiting the 
known scaling behavior of the susceptibility in conjunction with various hy- 
perscaling relations. We begin by computing the anomalous dynamical scaling 
dimension t] which corrects z at order 1/A^. 



4-4- 1 The anomalous dimension f] 

It is known that the uj = susceptibility should scale with momentum like 
G~^{k,0) ~ k^ where the bare dynamical critical exponent z = 2 will be 
corrected by the critical exponent r] as z = 2 — rj. Therefore, we can write 

G-\k,0)r~.k^~^ ,:^k^fl + ri\nj) (4.27) 

where A is a large momentum cutoff. From Eq. (14.111) at r = r^, T = and 
uj = we have 



^ ' ^ N J 27r'^'i 71 no(l,e. 



0) 



1 1 



;i + l/g)2 + e l + e_ 

(4.28) 
where we have used Eq. (ID.Sp and all variables of integration are dimensionless. 
Expanding the integrand for large q and identifying the logarithmic prefactor 
leads to 

1 7^ 3-e 

ae- 



' TT^iVy no(l,6,0)(e+l)3 

0.13106 ,, , 

^ ^j^- (4.29) 

Knowing the value of rj will be particularly useful because it will fix the cutoff 
dependence of the quantum critical conductivity at order 1/A^, since z = 2 — r], 
and we expect cr(T) ~ T^^/^ . Thus, if (y{T) = A/ \/T where A is a constant at 
A^ = oo, then at order 1/A^ we should have a = {A/VT)[1 + (r//2) ln(A/VT)]. 



4-4-2 The correlation length exponent v 

Calculating the 1/A^ correction to the correlation length exponent v is unfor- 
tunately not so simple, but we begin by examining the behavior of the inverse 
susceptibility at T = and /c = a; = as one tunes the coupling constant g 



56 



towards Qc 



6-1(0,0) -((?-(?, 



(4.30) 



which defines the susceptibihty exponent 7. At A^ = oo we know 7 = 2, and 
thus for finite A^ let us parameterize 7 = 2(1 — a), which can be related to v 
via the scaling relation 7 = (2 — r^)//. To this end, we define r^ via 



1 


1 f dk f duo 
g ~ J 2n J 27r 

_ v^ 


/ 1 
u' + kl 


1 


9c 


A;2+ uj\ + rg 



IT 



where we have exploited the fact that 



(4.31) 



27rno(0,0,rs 



(4.32) 



Thus, from Eq. (I4.3ip we have r^ ~ io ^ 9c)'^y and upon comparison with 
Eq. (H:30|) we find 



G-\0,0) r^ (g - g. 



\2(l-a) 



~ Tg ° ~ Tg 1 + a In ^^ 1 . 



A2^ 



(4.33) 



So again we can extract a critical exponent by determining the prefactor of a 
logarithmic divergence of G^^. At this stage it will useful to quote the following 
two results (with details given in Appendix [D| 



Ilo{k,uj,r) 
1 



2n\k\ 



k^ + \uj\ 



arcsm 



Y(P + |cj|)2 + 4Pr^ 



+ arcsin 



k^-U 



(P + |cj|)2 + 4Pr^ 



+ 



1 



An^k^ + 2\uj\ +4r 
-In 



In 



2Jr + IcuU/P + 2|cj| + 4r + P + 3|u;| + 4r 



2^r + Ic^l^P + 2|u;| + 4r - P - 3|c<;| - 4r 
2v/r./F + 2|u;| + 4r + P + |u;| + 4r 



2A/r + IcuK/P + 2\uj\ +Ar - k^ - \uj\ - 4r| 



(4.34) 



and 



no(A;,w,r) = 



(9no(/i;,u;,r) 



-2 



dr 
dq f de 



1 



2ti J 27T (g2 + |e| + r)2[(A; + q^ + |u; + e| + r 



■. (4.35) 
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Combining Eqs. (14.51) with (14.141) and (I4.3ip r can be written in terms of Vg 
which when used in Eq. (14.111) leads to the result for the inverse susceptibility 



2^v^ f dk f duj f dq f de 1 



G-i 0, 0) = r„ + F(r„) — ^^ / — / — / — / 



1 



A^ J 2tt J 2n J 27r J 27r no(g, e, 0) (P + |cj|)2 
1 1 



g2 + e (^k + qy + \uj + e\ 



(4.36) 



where we have been able to replace r with Vg in any term that is already of 
order 1/A^ and 



F(r„ 



1 r dk r du 



N J 27T J 27TUo{k,uj,rg 



1 U'Q{k,uj,rg 

k^ + \uj\+rg 2no(0,0,rg 



(4.37) 



A useful check is to note that G ^(0, 0) = above for Vg = 0. The next step 
is to investigate the small Vg behavior of F{rg). For this, let us first examine 



IIq as Tc, — i> 0, we find 



no(fc,w,rg) 
where 



1 ^ \uj\ 



1^ (\u:\ 



7r(P + |cj|)^ \k\^ \k^ J k^ "\k'^ 



^g + 



(4.38) 



$i(x) 



dy 



l/[4y2_6y+l + (6y-5) 



x\ 



47r(l + a;)y {2y - If \ ^/r^{y + xf/'^ 

4y^ - lOx^ + ll|/2 - 6y + 1 - (6y3 - 9y2 + Qy - i)x ' 



[y(l-y + x)]3/2 



$2(a;) = - 



2 3 + a; 



vr 1 + x 



(4.39) 
(4.40) 



and we have used the fact that |e| <C a; over the regime important for small 
Tg. From this expansion we can also determine the small Vg expansion of IIo 



Ilo{k,uj,rg 



Mk,u,0)- 



2./r: 






(4.41) 



58 



and finally that of F{rg 



T^Jrl f dk r duo <l>i(|a;|/A;^) r^ r dk f du 



^'^^^ N J 271 J 2n\k\mo{k,u;,0)^ N J 2tt J 27r 



2$i(|^IA^ 



+ 



7r$2(|cu|/P 



1 



|A;|3(A;2 + |cj|)n2(fc,w,0) fc4no(fc,w,0) (P + |tu|)2no(A;, w, 0) 

(4.42) 



Now, comparing this result with Eq. (14.331) and (14.361) the second term, which 
is linear in r^ defines a by 



F{rg) = --- + arg\n[—\ + 



(4.43) 



where a can be related to the correlation length exponent by the scaling 
relation 



7 = 2(1- a) = y{2-r]). 



Using Eq. (ID.SP to define 



(4.44) 



\UJ\ 



$ol^l =|A;|no(A;,a;,0) 



(4.45) 



and from Eq. (I4.42p a is given by 



a 



2Tr^N 
0.455 



duo 



2<l>i(u;) 7r<l>2(u;) 



(uo + l)^U^) ' ^q(uo 



[l + uoy^oH 



N 



(4.46) 



The value of u can finally be determined using Eq. (I4.29P as 



^ 0.389 
~ 1 . 



(4.47) 



The values found in this section for z = 2 — rj (Eq. (I4.29P ) and v (Eq. (14.471) ) 
corresponding to a A^ component complex field are fully consistent with previ- 
ous calculations including an e expansion near 2 dimensions (J65l : l50l ) (Eq. (I2.10p 
and fl2.1ip') and via Monte Carlo simulations where z = 1.97(3) and u = 
0.689(6) (IBJ. 
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4-5 Quantum transport at finite N 



We now endeavor to compute the dc values of the thermal and electrical 
conductivity in a 1/A^ expansion in the quantum critical regime. Transport 
is again calculated via the Kubo formula, and the initial steps are identical 
to those presented in Section [3.5.11 for the derivation of Eq. (13.751) . However, 
unlike the case where the number of components of our order parameter field 
was infinite, we now have the modified propagator of Eq. fl4.2ip and the single 
polarization bubble diagrams will be corrected by additional loops giving rise 
to corrections of order 1/A^. 

It will turn out that although the individual values of the thermal (k) and 
electrical (a) conductivities are not by themselves universal to order 1/A^, 
their ratio is a universal number, solely as a result of the appearance of an 
anomalous dimension that alters the critical dynamic scaling. 



4-5.1 Diagrammatic expansion 

We ask the reader to recall Eqs. (13.671) to (I3.75P and begin by writing down 
the expression for the transport coefficients (p = for electrical conductivity 
and p = 2 for thermal conductivity) obtained from the Kubo formula in terms 
of a polarization function at external imaginary frequency iun 

4g*2-p 

GpiiuJn) = -—ICpiiUn). (4.48) 

The current-current correlation function is evaluated with respect to Sr, Eq. (14.91) . 
and the result contains both a dia and paramagnetic part 

T r Hk r 

ICpizOO^) =- Y. TT (en +^n/2)^ [mk.en)?) 

- A;'(^:(fc, e„)^,(fc, e„ + u;„)^:(fc, e„ + ^„)^,(A;, e„)) 

As we are only interested in the real dc thermal and electric transport coeffi- 
cients, an examination of Eqs. (I4.48P and (14.490 indicates that after analytic 
continuation to real frequencies we will need the imaginary part of /Cp, and can 
thus focus only on the paramagnetic contribution to Eq. (14.490 corresponding 
to the four-point correlation function. This term arises as a result of quadratic 
insertions of the scalar or thermal potentials Aj. The resulting paramagnetic 
polarization function has the diagrammatic expansion to order 1/A^ given by 




ICl^-'^iuj^) = 6 h + 6 \ h (4.50) 
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where 



^ (4.52) 



and an open circle indicates a factor of /c(e„+ct;„/2)P/^ where p = for electrical 
transport and p = 2 for thermal transport and a closed dot represents the 
interaction i. The dashed line is the fluctuation propagator Ut/N and i?i is 
the finite shift in the critical point to order 1/N given by Eq. fl4.24p . The 
self-energy is defined in Eq. fl4.22p such that S(0, 0) = 0. To identify the role 
of various 1/A^ corrections to transport it will be useful to present the full 
integral form 

dk 



la^^'^iu^) =TY.j ^k\en + uj2fG^{k, en)G,{k, e„ + 

l^k\en + Un/2YGl{k, en)Go{k, e„ + ^„) 

/dh" 
— k\en + uj2YGl{k, er.)G,{k, e„ + uj^)E{k, e„) 



1t2 V [^ fdkk{k + q){e^ + uj2fl\en + Vt^ + u:j2Y/^ 



N ^^J 2nJ 27r Uriq^n^^R) 

X G^ik, en)GQ{k, e„ + u;„)G'o(^ + "?' e„ + ^n)GQ{k + g, e„ + n„ + cj„). 

(4.53) 

The first term is just the paramagnetic contribution in the large- iV polarization 
function previously defined in Eq. (I3.76p . The second term arises from the 1/A^ 
correction to the mass R, and the final two terms from the self-energy and 
vertex corrections respectively. 



4-5.2 Frequency summations 

The Matsubara summations in the various terms of Eq. (I4.53P can be per- 
formed by solving integrals in the complex plane with repeated use of the 



basic identity (l9ll ) 



T E -^(^^n) = IJ^, ^oth (^) [F{e + ^r/) - F{e - zr^)] (4.54) 

where we note that if F{ien) = !F{\en\) then after analytic continuation F(e ± 
irf) = ^(=pze). The full details on the derivation of various summation formulae 
used in this section are given in an Appendix [El The general approach will 
be as follows: use the relevant summation formula to obtain an expression for 
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each term in Eq. (14.530 analytically continued to real frequencies. Since we 
are only interested in dc transport, an examination of Eq. fl4.48p tells us that 
we will require the imaginary part of the term that is linear in the external 
frequency, a;. Thus by Taylor expanding our analytically continued result, we 
can extract the relevant transport coefficients. We will examine each term 
separately. 

For the first term in Eq. fl4.53p we could just as easily perform the Matsub- 
ara sum using the spectral representation of the bare Green function, which 
was done in Section 13.5.11 and led to Eq. fl3.8ip . This allows for a test and 
benchmark of the contour integration approach. We need to evaluate: 



^ {k'^ + R+ |e„|)(P + i? + |e„ + ujn\) 

= T^ J'2,p(«en,«(en + ^n)) (4.55) 



where we have suppressed the momentum dependence of T2,p for compactness. 
Using Eq. (IE. 30 we find 

1- Im/2,p(^ + ^r/) _ 1 f de e^+P 1 ..^^. 

o^^o u 2tJ 2n sinh2(e/2T) [(F + R)^ + e^]^ ' ^ -^^^ 

which does indeed agree with our previous result, Eq. (13.8 ID . Substituting into 
Eq. (I4.48P and defining 

Qp = lim Qp{iu + ir]) (4.57) 

= 0^=°° + ^f + Gp + Ql (4.58) 

where we have broken the total dc transport into a sum of four contributions 
coming from the four types of terms in Eq. (I4.53p . Because we have ignored the 
diamagnetic part of the polarization function, IC^'^^^ is purely imaginary and 
thus after analytic continuation Qp is a real number. The N = oo contribution 
is 

4e*^"P r de e^+P r dk k^ 



^=^ ^ 4e- - r de 6-^ r 



212 



2rp+i J 27r sinh'(e/2r) J 2tt [{k^ + R)^ + e 

1 f 0.217997- ■• e*^ ■,p = 
^{ . (4.59) 

V T 0.24592 • • • ; p = 2 



Due to the finite shift in the critical point, coming from Ri ~ 0(1/ N), we 
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need to evaluate a correction of the form 



h,pi'^^n) =TY^ 



:^ik^ + R+\en\yik^ + R+\en + uJn\) 

= T^^2,p{^('n,iien + UJn)) (4.60) 

en 

however, upon examination of Eq. fl4.55p it is clear that in the dc limit, this 
can be evaluated by taking a derivative of Eq. fl4.59p with respect to R. 



^'' - -im [(-^«')5'' 



N=oo 



1 I 0.062251- ■■ e*2 ;p = 

' (4.61) 



VTN [ 0.026867- ■■ ;p = 2 

where we have used the previously calculated values of R/T = 0.6248 and 
Ri/T = 0.1069/Ar. 

Examining the third term in Eq. (14.531) we now have to perform a dual Mat- 
subara sum over a function with four separate frequency arguments 

(4.62) 
= T^ Y^ J='4^p{ien,i^n,i{en + ^n),i{en + ujn))- (4.63) 

Using Eq. (1E.8P we can write: 

,. Iml^piu + ir]) 
lim = 

1 rdn rde t^e^+^k' + ^)csch2 (^) Re [Uriq, n, R)]-' 



(4.64) 



T J 2tI J 27r [(P + i?2)2 + g2]3 

(e + fi)coth(^) ficoth^2^ 



[(A; + g)2 + i?]2 + (e + fi)2 (^2 + ^)2 + ^2 
which leads to the self-energy corrections to the dc conductivities 
s _ 86*^-^2^ rdq_ rdn rdk r de e^^^k^k^ + R)csch^ {^) 

^P - fqXP+l J 2^ J 2^ J 2^ J 2^ [(P + R2y + g2]3 

(e + 1]) coth (^) f]coth(§) 



X Re[nT(g,fi,i?)]"^ 



[(A; + g)2 + /?]2 + (e + ny (g2 + /?)2 + ^2 j • 

(4.65) 
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The final term in Eq. fl4.53p has five separate frequency arguments 

1 



[(A; + g)2 + i? + |e„ + ^n\][{k + qY + R+ |e„ + n„ + cj„|] 

(4.66) 



Using Eq. (1E.14P we can write: 



1 /• df] /• rfe ^"[e(e + ^)\^^^l''c^c\i^ [^) csch^ (^) Re [n^lg, fi, i?)]"^ 



2Ty 27ry 27r [(F + i?2)2 + e2]2{p + ^)2 + ^]2 + (^ + ^)2}2 

X /(A;2 + i?)(e + fi)sinh('^) + [(A; + g)2 + /?]esinh(^^^)| (4.67) 

giving the vertex contribution to the o; = transport coefficients 

^__4e*2-^iP tdq r dQ r dk r de k{k + q)[e{e + Q)]^+p/^ 
^P - ~ jvTP+i J 2^7 2^7 2^7 2^ [(P + i?2)2 + e2]2 

csch^ (2^) csch^ (^) Re [nT(g, fi, R)]'^ 
"" {P + g)2 + i?]2 + (e + fi)2}2 

X |(A;2 + i?)(e + fi) sinh (^^) + [(A; + qf + R\esmh (^^) } • 

(4.68) 



4-5.3 Numerical evaluation 

The 1/A^ corrections to thermoelectric transport coming from the self-energy 
and vertex corrections are written in Eqs. (I4.65P and (14.681) as two four di- 
mensional integrals that cannot be evaluated analytically. Before we attempt 
to compute them numerically, we first present a simple argument concerning 
their expected ultra-violet behavior. From scaling we understand 



^p - ^ (4-69) 
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where in Section 
write 



we found that z = 2 — t] with rj ~ 0{1/N). Thus we can 



1 



rl/(2-r,) 



^'--{^^^-hjf 



(4.70) 



where Cp are universal constants and A is a non-universal ultra violet cutoff. 
Immediately we see that to order 1/A^, the ratio of the thermal to electrical 
conductivity divided by temperature — the Wiedemann-Franz ratio — will 
be independent of any cutoff as A ^- oo: 



W = 



Go 



gN=oc 



gN=oo 



'0 



1 + ^ + 



T 



A 



(4.71) 



This is an important equation that guarantees the universality of our final 
result, and will allow us to test the accuracy of our numerical integration 
procedure. 

We begin by combining the expressions for the self-energy and vertex correc- 
tions such that Re [11^(5, n, i?)]~^ (the most costly function to compute, as 
described in Appendix |D]) is in the outermost integral. 



g"" + g' 



dk 



dq 

duj 

2^ 



de 

2^^ 



Re 



UT{q,n,R) 



(4.72) 



Y''{k,e,q,n) + Y^^{k,e,q,n) 



where 



V^(fc,e,g,n) 



2k\k^ + R)e^+Pcsch^ [^) 

[(P + Ry + e2]3 
+ fi)coth(^) ficoth(|.) 



Yl'{k,e,q,n) 



[{k + qY + Rf + (e + f])2 (g2 + R)2 + ^2 

k{k + q)[e{e + 9)Y+Pl\^di^ (0) csch^ ( 



(4.73) 



2T 



p2 + ^)2 + ,2]2|p + qY + ^]2 + (, + ^yy 



{k^ + R)ie + n) sinh {-) + [{k + qf + R]esmh 



T 



(4.74) 



Performing the outermost integral numerically using an adaptive routine we 
arrive at the final results shown in Fig. [16] where a tilde indicates that a 
quantity has been multiplied by a factor of A^vT/e*^~^. After fitting to the 
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Fig. 16. The 1/N corrections to the rescaled thermoelectric transport coefficients 
coming from self-energy and vertex corrections {p = 0, bottom curve for a and 
p = 2, top curve for k/T) plotted as a function of a dimensionless external ultra 
violet momentum cutoff A/yT- The solid lines are fits to the expected divergent 
behavior, from which the non-divergent corrections as A ^ oo can be extracted. 

expected divergent form in Eq. (14. 70 p . we find (as A ^ oo) 



^ r e*2 0.118 ^^^r/ A 



gf + Q. 



Vt N 

1 0.151 



-,N=ooV 



In 



A 



T 



(4.75) 
(4.76) 



and combining with the previous results of Eq. (14.591) and (14.611) we have the 
dc thermoelectric transport coefficients to order 1/A^ 



.*2 



a 

K 

T 



T 



0.218 + 



0.0561 0.0142, A 



A^ 



+ 



0.246+^ + 



A^ 
0.0161 



In 



T, 



/T \ ' N N 

which both explicitly depend on A as expected. 



In 



A 



T, 



(4.77) 
(4.78) 



4-6 Wiedemann- Franz ratio in the quantum critical regime 



We have evaluated the full fluctuation corrections to thermoelectric transport 
up to order 1/A^ (Eqs. (14.771) and (I4.78P ) coming from the direct contribu- 
tions of Cooper pairs due to the proximate superconducting state. Initially 
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Fig. 17. The 1/A^ corrections to the Wiedemann-Franz ratio plotted as a function of 
the inverse external rescaled ultra violet momentum cutoff in the quantum critical 
regime. The solid line is a fit to a second order polynomial and the individual diver- 
gences of the electrical and thermal conductivity are exactly canceled as A/yT — > oo 
giving the universal correction 6W = 0.0376/A^. 



dismayed by their cutoff dependence, we now recall the previous argument 
that led to Eq. (14.711) . We found from scaling that the required T~^'^ tem- 
perature dependence of k/T and a implied that when dividing them to form 
the Wiedemann-Franz ratio, all divergent A-dependence must exactly cancel. 
This exact cancellation is seen in Fig. flT7|) where we plot the total correction 
to the WF ratio, 5W as a function of the inverse of the rescaled dimensionless 
cutoff. As A — > oo, 5W approaches a constant. Extracting the infinite cutoff 
result via a polynomial fit we find (after inserting the proper power of the 
Boltzmann constant) 



W 



of 



0.0376\ fke^ 



(4.79) 



Therefore, the WF law is indeed obeyed, (i.e. is temperature independent) 
indicating the presence of only fully elastic scattering and is independent of 
any microscopic constants. The term proportional to 1/A^ is quite small, and 
for the physical case, A^ = 1 it corresponds to a correction on the order of ten 
percent. 
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5 Conclusions 



This paper has been concerned with a topic that could be mistakingly con- 
fused with one of limited scope, the pairbreaking quantum phase transition 
between a superconductor and a metal in an ultra-narrow wire as modeled by 
a continuum quantum field theory. Instead, we have discovered a remarkably 
rich phase diagram full of interesting phases and crossovers. 

Experimental motivations for a theoretical analysis of this transition exist 
in the form of transport experiments on metallic nanowires, formed through 
molecular templating by sputtering material on top of a long rigid "bridge" 
or "backbone" molecule lying over a trench (l2|). In this way, wires with diam- 
eters of less than 10 nm can be fabricated; a giant step towards reaching the 
limit where the length scale characterizing quantum fluctuations approaches 
the finite radius of the wire. In an applied current and at fixed temperatures, 
below the bulk superconducting transition temperature for the wires com- 
posite material, a given wire can display either metallic or superconducting 
behavior depending on its radius, with the general trend that thinner wires 
are less superconducting. In addition, for a particular wire which does exhibit 
electrical transport without resistance, superconductivity can be destroyed by 
turning on a suitably strong magnetic field oriented along its parallel axis. In 
both of these cases, it is some non-thermal parameter which suppresses the 
Cooper pairing instability and tunes between the superconducting and normal 
metallic state at zero temperature, providing an excellent manifestation of a 
quantum phase transition. 

The description of the transition that we have adopted in this paper is a critical 
theory of strongly repulsive, fluctuating Cooper pairs, written in terms of a 
complex order parameter, overdamped by its coupling to a bath of unpaired 
fermions resulting from the presence of some type of pairbreaking interaction. 
The existence of the bath, imagined as a large number of unpaired electrons 
residing in the transverse conduction channels of the wire, leads to a long 
range interaction in imaginary time providing Ohmic dissipation in the form 
of a non-analytic |u;„| term in the effective action. The presence of such an 
anisotropic relationship between space and time coupled with a continuous 
symmetry order parameter fixes the dynamical critical exponent at z = 2, 
and the resulting upper critical dimension is (iucD = 2. 

The thinness of the experimentally investigated wires provided us with a useful 
theoretical framework in the form of the quasi-one dimensional limit, where 
the radius of the wire R is on the order of, or much smaller than the su- 
perconducting coherence length ^ at low temperatures. The length scale ^ 
measures the average separation between the electrons in a Cooper pair, and 
ioT R < ^ the composite object begins to feel the boundary. If the wire is suf- 



68 



ficiently long, the paired states can be described in terms of a quantum field 
theory in one space and one imaginary time dimension. In 1 + 1 dimensions, 
we found ourselves below ducr,, with the repulsive self- interactions between 
Cooper pairs being strongly relevant. As a result, any perturbative or mean 
field approaches are unable to provide a complete and accurate picture of the 
physical phenomena. 

To deal with strong interactions, we have employed a variety of field theoretic 
and numerical techniques suitable in different regions of the pairbreaking phase 
diagram in conjunction with a careful scaling analysis. In the quantum critical 
regime, at temperatures limited to those where hj -C k^T, dynamical observ- 
ables were computed from the real time Langevin dynamics of an effective 
classical theory renormalized by quantum fluctuations in the limit where the 
number of complex order parameter components {N = 1 in the physical case) 
was equal to infinity. Near the finite temperature classical phase boundary, 
the effective theory was used to derive a modified version of the thermally 
activated phase slip theory of Langer, Ambegaokar, McCumber and Halperin 



( llll : Il2l ) which does not explicitly depend on the number of transverse chan- 
nels of the wire. At low temperatures, the full quantum theory was analyzed 
in a systematic expansion to first order in 1/A^, with extensive detail provided 
on the calculation of critical exponents and various physical quantities near 

z = 2. 

The combination of all these approaches allowed us to determine the full form 
of the zero frequency (dc) electrical and thermal conductivities as a function 
of temperature and the pairbreaking parameter which drives the transition. 
Our first experimentally testable result is a complete crossover phase diagram 
for the quantum superconductor-metal transition (SMT). We predict that 
upon fixing the source of pairbreaking (either magnetic field or wire radius) 
at a value near criticality, as function of decreasing temperature, the electrical 
conductivity should first increase as 1/VT due to pairing fluctuations near 
the quantum critical point, then change to decreasing as T^ once the low tem- 
perature metallic phase has been reached. The T~^/^ dependence of the dc 
electrical conductivity appears as a natural consequence of scaling near criti- 
cality and was not found in previous non-interacting microscopic perturbative 
calculations (44l). There may already be qualitative experimental evidence for 



transport that is non-monotonic in temperature (J82l : l33l ) and a further analysis 



of the reported resistance data along these lines seems apropos. 

The second prediction is that in the quantum critical regime at finite temper- 
atures, defined by a pairbreaking strength that is close to the one that would 
destroy superconducting order at zero temperature, the ratio of dc thermal (k) 
to electrical (a) conductivity divided by temperature (the Wiedemann-Franz 
ratio) should be a temperature independent constant measuring deviations 
from the Lorenz number for a normal metal /q = 7r^/3(A;B/e)^. We have com- 
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puted the exact value of this correction in a 1/A^ expansion and found 

^^ fr.^^^ 0.0376\ fkB^'^ 



».._. ^0.282 + ^^j f . (5.1) 



Although electrical transport measurements already exist on a variety of metal- 
lic nanowires, their thermal contact with both the substrate as well as the two 
dimensional leads remains a challenging problem in the suspended molecular 
geometry, and would inhibit experimental access to this prediction. However, 
various intriguing and promising directions are currently being explored in- 



cluding the use of a scanning tunneling microscope (J92l ) 



Further avenues for theoretical progress still remain in superconducting sys- 
tems in confined geometries, including a full understanding of the superfluid 
density in the low temperature ordered phase which has not been attempted 
here. Such a description would require proper inclusion of the pairing inter- 
action as well as Coulomb repulsion in the presence of disorder, leading to 
a plasmon mode describing the strongly fluctuating phase of the supercon- 
ducting order parameter. In addition, the presence of a recently identified 
infinite randomness fixed point at the disordered superconductor-metal tran- 



sition (J72l : l93l ) opens up a plethora of questions surrounding the nature of 
hydrodynamic transport in the presence of unconventional activated dynamic 
scaling and a detailed investigation by the authors is currently underway. 
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A Computation of /(f2) 



This appendix provides details on the evaluation of the scaling dimension 
of p characterizing the strength of particle-hole symmetry breaking in the 
perturbing action Sp. The combination of the two-loop diagrams in Eq. (12.141) 
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evaluated at zero external frequency and momentum lead to 



^(")— Vi^/i^/rl^/ 



{27iY J {27tY 



(A.i; 



(P + |a;i|)2(g2 + |^2|)[(fc + qy + |^i + ^2 + ^1] ' 



The two momentum integrals can be evaluated by employing Feynman pa- 
rameters 



I(Q) = -i2pu 



(4vr) 



1 1 
dx / dy 



,TiA-d) } } fdu, fduj2 (i-yy-"/2 







2n J 27r [1 - y{x^ -x + 1)^/^ 

(^ 

[(1 - y)\uji\ + y{l - x)\uj2\ + xy\uji + uj2 + n\]^~'^' 

(A.2) 



where T{x) is the gamma function. The double frequency integral, can be done 
by determining the sign of the various absolute values in the relevant seven 
regions of the uJi — uj^ plane delineated by changes in signs ojx , 002 and ujx + 002 
leading to the following useful but complicated looking result 



dijj\ r duj2 



UA,B,C..)^I'^I 



iVi 



27r {A\uji\ + B\uj2\ + C\uJi + uj2 + fi|)'" 



ABCT{(r - 3)1]3- 
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a 
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A2(A2-C2)(A2-52) 
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2^2 



(52-C2)(A2-52)2 

(A.3) 



In (i = 2 — e and a = 2 + e, J,^ has a pole at e = with a residue that can be 
read off from 



J^(A,5,C,2 + e) = -^ 



5C(2A + 5 + C) 



(A + 5)2(A + C)2(5 + C) 

Using Eq. (lA.4p in Eq. (1A.2P we find 



- + 0(1). (A.4) 
e 



xy{\-x){\ -y){2-y) 



i^u p r r 

^^^^^^^1 "^""l '^^ir^y(x2-X+l)][l-l/ + Xy2(l-x)]2 



(A.5) 
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B Connection to microscopic BCS theory 



In order to motivate the experimental relevance of the effective action Sa, 
the microscopic values of the renormalized pairbreaking frequency a, the bare 
diffusion constant D, dissipation strength 7 and quartic coupling u can be 
determined in both the clean and dirty limits. We begin with the connection 
of the pairbreaking frequency to various experimentally relevant geometries, 
then move on to the relationship between the theory presented here and a 
time dependent Ginzburg-Landau (TDGL) theory for a conventional super- 
conductor. 



B.l Pair-breaking in quasi-one dimensional wires 



As mentioned in the introduction, there are various origins of pairbreaking 
perturbations relevant to experiments on ultra-narrow wires. The most theo- 
retically appealing consists of magnetic impurities localized on the surface of 
a metallic wire leading to an inhomogeneous BCS coupling (Fig. [2]). In this 
case, the microscopic value of a is not known exactly, but it can be related 
to the inverse of the spin-flip scatting time. However, there are a number of 
well-defined experimental geometries where the actual value of a can be com- 
puted in terms of the physical properties of the system. One case which is of 
interest here, is a narrow metallic wire in a parallel magnetic field. 



In the dirty limit. Shah and Lopatin (J44l ) have computed the precise form of a 



through the use of the Usadel equation formalism. They find that for a narrow 
diffusive wire with radius R smaller than both the superconducting coherence 
length and the magnetic penetration depth placed in a parallel magnetic field 



a^ 



D /eHR\^ 



) (B.l) 



2 V c 
where D is the diffusion constant and c the speed of light. 

B.2 Microscopic parameters in the clean and dirty limits 



The microscopic values of D ,7 and u can be found through an analysis of 
the time dependent Ginzburg-Landau theory studied by Tucker and Halperin 



(J94J ). There, the three dimensional equation of motion for the Cooper pair 



operator \E'(a;,t) in real time is given by 

d_ 

"dt 



d r 1 

%— ^(a;,t) = - a + b\<i/{x,t)\'^ + 6{-iVf ^(a;,t). (B.2) 
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Rescaling, to ensure that the coefficient of the time derivative term is unity 
and performing an integral over the cross-sectional area of the wire to move 
to the quasi-one dimensional case of interest {'^{x,y,z,t) ~ \l/(x, t)), we read 
off the value of the coupling constants to be 



D 



u 



_5_ 



AhW 



(B.3) 
(B.4) 



where A is the cross-sectional area of the wire. Appendix A of Ref. (J94j ) gives 
the microscopic values of 5, b and 7 as 



2m 



h 

7 



h^ 



2me(0)nx(0.882eo//)^ 



rn^ 



lQme{Q)kB%o 



(B.5) 
(B.6) 

(B.7) 



where ^(T) is the Ginzburg-Landau coherence length, ,^0 the BCS coherence 
length, £ the mean free path and x{p) the Gor'kov function defined by 



x(p) = E 



n=0 



{2n+lf{2n + l+p) 



E 

.n=0 



(2n+l)3 



(B.8) 



The critical temperature (Tc) and density of conduction electrons in the normal 

state in) are known to be 



eo 



ksTcQ 
n 






37r2 



(B.9) 
(B.IO) 



respectively, while, the zero temperature coherence length and relevant Gor'kov 
function depend on whether we are in the clean or dirty limit 



m 



x(0.882eo//) 



0.74^0 ; 6 < 
0.85v^ ; ^0 » 

1.33£/eo ; eo > ^ 



(B.ll) 



(B.12) 



We have now gathered all the required information to compute the actual 
microscopic values of our model parameters in the dirty (^0 ^ ^) and clean 
(^0 ^ ^) limits differentiated by the subscripts d for dirty and c for clean. 
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Using the above relations, we find that for the dirty limit 

1 



Dd 
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D = 

1.5 


r^' 


Id 


r^ 


kpi 




Ud 


r-^ 


2.9- 


vf 

■ A T 



hN, 



(B.13) 
(B.14) 
(B.15) 



where the number of transverse conduction channels is assumed to be large, 
and is given by 

N^ = ^. (B.16) 

Similarly, in the clean limit 

Dc = \vF^o (B.17) 

2 n 

7c ^ TTT- (B.18) 

u, = Ud-2.9-£- (B.19) 

where we note that the bare value of the quartic coupling is identical in both 
limits. The value of these parameters clearly depends on the particular nor- 
malization scheme chosen for the order parameter, but the final results for all 
physically measurable quantities, such as the conductivity, will obviously be 
normalization independent. 



C Ginzburg-Landau theory near Tc 



The insertion of a self-consistent infrared cutoff allowed for the derivation of an 
effective potential at zero temperature within the ordered phase (Eq. fl3.52p ). 
At finite temperatures, near T^ Eq. fl3.37p can be expanded in powers of |crp 
in order to derive an effective Ginzburg-Landau theory for the superconduct- 
ing phase with coefficients renormalized by quantum fluctuations. This is ac- 
complished by first computing the critical temperature where (a) becomes 
non-zero. Suitable derivatives are then taken to determine the values of the 
quadratic (ao) ^t^ quartic {(3) coefficients in the expansion 

Vgl = Vo + «o(T - TJlap + ^f3\a\^ + ■■■ . (C.l) 
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Fig. C.l. The rescaled critical temperature found from the solution of Eq. ()C.3P as 
a function of /(e) given in Eq. ()3.49p which is shown in the inset. 

C. 1 Evaluation of the critical temperature T^ 



We begin by considering the saddle point equation (Eq. (13.331) ) in the presence 
of the symmetric cutoff. A similar procedure that led to Eq. (13.361) can be used 
here, giving 



I |2 ^I'^l ^9^ 
VD 7r2 



2 - In I ^^^ 1 + V^ ' 



2txT 



9 



vr^ 



dk 
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Dk'^ + r 



ij\l + 



27cT J 

DP + r" 
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In 



'DP' 
2-nT 



(C.2) 



Returning to Eq. (13.461) and noting that at T = Tc, fi\o'\^ = 0) = 0, we can 
derive an equation for T^ from Eq. ( JC.21) 



^Wr, 7r2/(e) 72^ 



2 -In 



2X2 



{l + efS 
2vrTJ2(e) 



+ ip 



2X2 



2vrTJ2(e) 



TX^ 



dk 



(1+^)1^1 



2P 



-?/'(l + A;2)+lnA;^ 



(C.3) 



where fc is a dimensionless momentum, and the cutoff can be fixed by setting 
e and using Eq. (I3.49p . Solving numerically using a secant method gives the 
result seen in Fig. IC.lt where /(e) is plotted as an inset. The critical tempera- 
ture is found to be proportional to the square of the distance from criticality. 
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Tc oc S"^ and more specifically 



Tc = ciie)5' 



(C.4) 



where ci ~ 1.90 for /(e) = 1. One could either fix e at this value, or choose a 
value of e using a plot like Fig. IC.ll that reproduced the relationship between 
Tc and 6 measured in an experiment. 



C.2 Evaluation of the quadratic coefficient ao 



In order to evaluate ao in Eq. ( JC.ll) we again appeal to Eq. f l3.46p and note 
that 

ao{T-T,) = -r{\a\' = 0). (C.5) 

9 



Near Tc, r <^ 1 and a double expansion of Eq. (1C.2I) can be performed in r and 
the reduced temperature t = {T — Tc)/Tc leading to (after some considerable 
algebra) 
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where only linear terms in r and t have been retained, ijj^^^x) is the first 
poly gamma function and 



1 
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(C.7) 



Comparing Eq. (1C.6P with Eq. (JC.SP we find 
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and upon choosing /(e) = 1, C2 — 0.856 and ao — 0.509385/(7. 
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C.3 Evaluation of the quartic coefficient f3 



In order to determine the value of the quartic coefficient, we examine Eqs. (13.461) 
and (I3.53P at T = Tc leading to 
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(fViG 1 dr 
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T=T^ 



(C.9) 



Taking a derivative with respect to |crp of Eq. flC.2p we find 
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(CIO) 



and using the relation between A and ^(0), Eqs. (I3.48P and (I3.49P as well as 
Eq. ( IC.Qp . the quartic coefficient is given by 
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(C.ll) 



where Eq. flOTl) has been used. For /(e) = 1, /5 ~ 0A95VD\6\/g^ 



The results of Eqs. flC.3p . fIC.Sp and fIC.lip constitute all the required ingre- 
dients needed to construct the Ginzburg-Landau potential of Eq. ( JC.ip with 
the final result shown in Fig. [91 This appendix has presented numerical results 
for a fixed value of e that was determined by setting /(e) = 1. However, the 
/(e) dependence of all coefficients is relatively weak for /(e) > 1 which is well 
satisfied by a physical value of the cutoff. 



D The Fluctuation Propagator 



In this appendix, we will provide details on various results related to the 
evaluation of the fluctuation propagator at both zero and finite temperatures. 
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D.l Zero temperature 



At zero temperature, and coupling g the fluctuation propagator is given by 



dq f dt 
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The momentum and frequency integrals can be done by employing Feynman 
parameters to yield 
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Using the relation 



Re [atanh(2;)] = - In 
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we can write 
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and at the critical coupling where r = this simplifies to 
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D.2 Finite temperature 



A key step required for the evaluation of the shift in the critical point coming 
from 1/A^ corrections (Eq. (14.241) ) and the thermoelectric transport coefficients 
(Eq. ( J4.72P ) is the fast and accurate computation of IlT{k,Un,R) as well as 
the real and imaginary parts of its analytically continued value just above the 
real axis where iun -^ u + irj. 



D.2.1 Numerical evaluation 



Starting from Eq. (14.71) and performing the momentum integral we have 
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Let us first find the value of the finite temperature fiuctuation propagator at 
k = ujn = Q- Starting form Eq. (1D.6P above we can derive a simple result 
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where ({m,x) is the Hurwitz zeta function. In order to evaluate the sum in 
Eq. (1D.6P at finite frequencies and wavevectors, we explicitly sum the terms 
up to some large value of |e„| < 2nL where 27rL 3> Ic^^l. For the remaining 
terms, we perform a series expansion of the summand in powers of l/|e„| and 
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then use the asymptotic series 
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As discussed in Ref . (J90l ) , we must use the value of R given in Eq. fl4.17p for 
the resulting Urik, Un, R) to be well behaved at large k and Un- 



D.2.2 Re[nT(g,n,i?)]-i 

We now provide details on the use of the summation formulas described in 
Appendix [E] to evaluate the real and imaginary parts of the fluctuation prop- 
agator analytically continued to real frequencies. The benefits of this rather 
complicated derivation are manifest in the increased computational efficiency 
of having analytic expressions as opposed to resorting to a Kramers-Kronig 
relation. Again, we start from 



UT{q,nn,R) = TY.I 



dk 



1 



27r {k^ + R+ \en\)[{k + qy + R+ |e„ + fij] 



(D.9) 



and need to analytically continue to real frequencies ifi„ — > Q + irj. Thus, 
using Eq. (IE.3P we perform the Matsubara summation to give 
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where 
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Considering each of the four terms in Eq. (1D.10|) separately, we will have to 
perform an integral of the form 
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dk 
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(D.12) 
which was evaluated using Feynman parameters and in the particular case 
considered here Re a = Reb = R > 0. Using the relation 
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Im 
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q'^ + iVa- Kbb + V« - Kccf A{q,a,bXb,c,Cc) 



Cbb + CcC + Cc sgn(c) V(v/^^T62 + a)iV^^T^ 



+ CfeSgn(6)V(Va2 + 62-a)(Va2+^ + 



(D.15) 



where 



A{q,a,b,Cb,c,Cc 
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+ 



Cbb + CcC + Cc sgn(c)7(Va2 + 62 + a){^a^ + c^ - a) 



+ CfeSgn(6)V(\4^rr62-a)(v/^^T^ + 



(D.16) 



with a^b^c &M. and Cb.c = ±1- 



Therefore, using Eq. (1D.12P to Eq. (10.160 we can write the real and imaginary 
parts of the analytically continued fluctuation propagator as (suppressing all 



R dependence) 



RenT(g,l^±ir7) = 2/^1 ^^^^VXr) I'^/n(g,e + ^^,e + fi + ^^) 
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and 
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where (^,1/ = ±1- Such a formulation allows us to compute both the real and 
imaginary parts of IIt without having to resort to a Kramers- Kronig relation, 
leading to 
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UT{q,n,R) 



ReUT{q,n,R) 



[Re Uriq, fi, R^ + [Im Uriq, ^, R)Y 
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E Details on the Evaluation of Matsubara Sums 



This appendix provides details on the evaluation of multiple Matsubara sum- 
mations which appear in the evaluation of various current-current correlation 
functions coming from the Kubo formula. We begin with the basic identity 



CO 
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noting that if F{ien) = J?-'(|e„|), then after analytic continuation F(e ± irj) 
T{^it). By a similar application of contour integration we obtain 



h{i^n) =T^ J'(ze„,z(e„+a;„)) 



de , / e 
2 7 2;u^°*H^ 
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(E.2) 



and so 
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e + uj 
2T 



(E.3) 



where in the last expression we only denote the sign of the irj term, because 
the frequency arguments remain the same in all terms: 



F(±±) =F{e±ir],e + uj± ir]). 



(E.4) 



Rearranging the terms to preserve the order of the frequency arguments will 
allow us to pull out common factors in the numerator and lead to many 
simplifications. 

Any corrections coming from the presence of a finite self energy at order 1/A^ 
require that we perform a dual Matsubara summation over a function with 
four frequency arguments. Through a further generalization of the method of 
contour integration used to obtain Eq. (]E.2p we find 
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where 



T|(e, fi,-iu;, 



'n 



+ F{e-i' 



iuj- 



'n) 



it], Q + iri,e + iojn 
iri,Q — iri,Q + e — irj, e + iu. 
+ F{e — irj,VL — e + irj,VL — 



+ F(e 
+ F{e 



■ir],e + iuJn) 
- iuJn, Q + irjjQ + e 



and 



T4 (e, Q, iun] ir]) 



iujn,e + if]) 
■ iun, fl — e + iun, fi + iri,e + irj) 
+ F(e — ioJn, VL — irj,VL + e — iun, e — irj) 
+ F{e — iUn, fi — e + iUn, Vt — irj.e — irf) 



F{e + iri,Q — iri,Q + e + iri,e + iUn) 

+ F{e + iri,Q — e — iri,Q — iri,e + iuJn) 

+ F{e — iri,Q + iri,Q + e — iri,e + ioJn) 

+ F(e — irj.Vt — e + irj,VL + irj.e + ioJn) 

+ F(e — ioJn^ VL — irj,VL + t — iUn, e + irj) 

+ F{e — iUn, f2 — e + iUn, 0, — iri,e + irj) 

+ F{e — iUn, Q + iri,Q + e — iUn, e — ii]) 



so that 
h{iu + ir]) 



2TTi 



+ Fie 



de 

27ri 



tUr, 



coth I — - I coth ( — - 
\2TJ \2T 



+ coth 



n + e 

2T 



coth I — - 
V2T 



(E.6) 



fi — e + iUn, fi + ir],e — ir]) (E.7) 



+ coth — coth 

\2T) V 2T 



+ coth I _„ I coth 



2T 



2T 



F(+ + ++)-F(+-++) 

+ F{ +)_F(- + -+) 

F(+-++)-F(+--+) 

+ F(- + -+)-F(- + ++) 
F(- + -+)-F( +) 

+ n )-n- + —) 

F(- + ++)-F(- + -+ 



+ F(- + 



+ +- 



(E.8) 



85 



Multiple change of variable transformations have been performed to ensure 
that the function F in the terms above have the same arguments e, fi, e + fi, 
e + a;, and thus only the signs of the ir] terms have been denoted, 



F(± ± ±±) = F{e±iri,Vt±iri,t + Vt±iri,t + uj± ir]). 



(E.9) 



Finally, the vertex corrections are similar in that we still need to perform a 
dual Matsubara sum, but now the frequency arguments are more complicated, 
and there are five unique combinations 
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and 
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Again performing multiple variable shifts yields a much simpler expression 
where the frequency arguments of each term are the same. Suppressing the 



frequency dependence of /s = I^{uj + irf) we finally arrive at 
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where the arguments of F have been shifted tohe e, e + u, e + Q, e + Q + u, 

n, 

F{±, ±, ±, ±, ±) = F{e±iri, e+uo±ir], e + VL±ir], e + VL+uo±ir], VL±irf). (E.15) 
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